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Abstract. We consider the weakly asymmetric limit of simple exclusion process with drift to the 
left, starting from step Bernoulli initial data with p_ < p+ so that macroscopically one has a 
rarefaction fan. We study the fluctuations of the process observed along slopes in the fan, which 
are given by the Hopf-Cole solution of the Kardar-Parisi-Zhang (KPZ) equation, with appropriate 
initial data. For slopes strictly inside the fan, the initial data is a Dirac delta function and the 
one point distribution functions have been computed in [2J, [32]. At the edge of the rarefaction 
fan, the initial data is one sided Browman. We obtain a new family of crossover distributions 
giving the exact one point distributions of this process, which converge, as T / oo to those of 
the Airy Az^bm process. As an application, we prove moment and large deviation estimates for 
the equilibrium Hopf-Cole solution of KPZ. These bounds rely on the apparently new observation 
that the FKG inequality holds for the stochastic heat equation. Finally, via a Feynman-Kac path 
integral, the KPZ equation also governs the free energy of the continuum directed polymer, and 
thus our formula may also be interpreted in those terms. 



1. Introduction 

It is expected that a large class of one dimensional, asymmetric, stochastic, conservative interacting 
particle systems/growth models fall into the Kardar-Parisi-Zhang (KPZ) universality class. A 
manifestation of this is that the KPZ equation should appear as the limit of such systems in the 
weakly asymmetric limit. The weakly asymmetric limit means to observe the process on space 
scales of order e -1 and time scales of order e -2 , while simultaneously rescaling the asymmetry of 
the model so that it is of order e 1 / 2 . This sort of weak asymmetry zooms in on the critical transition 
point between the two universality classes associated with growth models - the KPZ class (positive 
asymmetry) and the Edwards Wilkinson (EW) class (symmetry) - and thus further confirms a 
mantra of statistical physics that at critical points one expects universal scaling limits. 



Bertini and Giacomin |12] obtained the first result for the weakly asymmetric simple exclusion 
process near equilibrium. This is extended to some situations farther from equilibrium in [2], 
directed random polymers in [T], and partial results are now available |21] for speed changed 
asymmetric exclusion. In this article we study the situation where asymmetry is to the left and 
the initial data has an increasing step, so that in the hydrodynamic limit one sees a rarefaction 
fan. We observe the process along a line x = vt within the fan and study the fluctuations. These 
converge to the KPZ equation with initial data depending on v. For v strictly inside the fan, the 
initial data is an appropriate scaling of a delta function, and the distribution of the fluctuations 
is known exactly [2], |32| . Our main interest in this article is the fluctuations at the edge of the 
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rarefaction fan. The scaling turns out to be a little different, but the fluctuations are still given 
by KPZ. Note that since the work of [12] it is understood that KPZ is only a formal equation for 
the fluctuation field and is rigorously defined as the logarithm of the stochastic heat equation. The 
edge fluctuations correspond to starting the stochastic heat equation with 

exp{-B(X)}l x>0 

where B(X) is a standard Brownian motion in X, with B(0) = 0. We will obtain an exact expression 
for the one point probability distribution of the resulting law - the edge crossover distribution - at 
any positive time. The main tool is the Tracy- Widom determinental formula for one sided Bernoulli 
data and therefore we are restricted to asymmetric exclusion. The resulting law is expected to be 
universal for fluctuations at the edge of the rarefaction fan for models in the KPZ class. 



1.1. Height function fluctuations at the edge of the rarefaction fan for ASEP. The 

asymmetric simple exclusion process (ASEP) with parameters p, q > (such that p + q = 1) is 
a continuous time Markov process on the discrete lattice Z with state space {0, 1} Z (the l's are 
thought of as particles and the O's as holes). The dynamics for this process are given as follows: 
Each particle has an independent exponential alarmclock which rings at rate one. When the alarm 
goes off the particle flips a coin and with probability p attempts to jump one site to the right and 
with probability q attempts to jump one site to the left. If there is a particle at the destination, 
the jump is suppressed and the alarm is reset (see [27] for a rigorous construction of this process). 
If q = l,p = this process is the totally asymmetric simple exclusion process (TASEP); if q > p it 
is the asymmetric simple exclusion process (ASEP); if q = p it is the symmetric simple exclusion 
process (SSEP). Finally, if we introduce a parameter into the model, we can let q — p go to zero 
with that parameter, and then this class of processes are known as the weakly asymmetric simple 
exclusion process (WASEP). It is the WASEP that is of central interest in this paper since it 
interpolates between the SSEP and ASEP, and is intimately connected with a stochastic partial 
differential equation known as the KPZ equation. We denote the asymmetry 

7 = q - p. 

We consider the family of initial conditions for these exclusion processes which are known of as 
two-sided Bernoulli and which are parametrized by densities p~,p+ £ [0,1]. At time zero, each 
site x > is occupied with probability p + and each site x < is occupied with probability p_ (all 
occupation random variables are independent). These initial conditions interpolate between the 
step initial condition (where p_ = and p + = 1) and the equilibrium or stationary initial condition 
(where p_ = p + = p). We will focus on anti-shock initial conditions where p_ < p + . 

Associated to an exclusion process are occupation variables r/(t, x) which equal 1 if there is a particle 
at position x at time t and otherwise. From these we define spin variables fj = 2rj — 1 which take 
values ±1 and define the height function for WASEP with asymmetry 7 = q — p by 



h-y(t, X) 



2N(t) + J2 Q<y < x v(t,y), x>0, 
2N(t), x = 0, 

{^(t)-E x<y <ov(t,y), x<o, 



where N(t) is equal to the net number of particles which crossed from the site 1 to the site in 
time t. Note that at time t = 0, /i 7 (0, x) is a two-sided simple random walk, with drift 2p + — 1 in 
the positive direction from the origin and drift 2p_ — 1 in the negative direction. 
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Figure 1. Limit profile for the height function of ASEP with two-sided Bernoulli 
initial conditions with /)_ = and p+ = 1/2. The rarefaction fan corresponds to 
velocities v < 0, and the edge to velocity v = 0. 



Proposition 1 (Hydrodynamic limit). Let p- < p+, 7 = e 1 / 2 and t = e 3 / 2 . Then, in probability, 

r 2p_(l-/9_) + (2p_-l)v forv<2p--l, 

< (l + v 2 )/2 forve[2p--l,2p+-l], (1) 

2p + (l-p + ) + (2p+-l)v forv>2p+-l. 



lim 



For 7 positive and not going to zero with e this result is well known [311 135] . We were not able to 
find a reference in the weakly asymmetric case. It is an easy consequence of the fluctuation results 
(i.e., Theorem 1 16j) which make up the main contribution of this paper (see however, Remark 1 17j) . 

The region v £ (2p- — 1, 2p + — 1) is the rarefaction fan, while v = 2p± — 1 is the edge of the fan. 
See Figure [1] for an illustration of this limit shape. 



For the purposes of this introduction let us set p- = and p + = 1/2 so that the right edge of the 
rarefaction fan is at velocity v = 0. Around this velocity one sees a transition from a curved limit 
shape for the height function to a flat limit shape. According to (pQ), the limit height is t/2. 

Definition 2. For m > 1, e > 0, T > and X 1} . . . ,X m G R set 

t = e- 3 / 2 T, x k = 2 1 'H 2 ^X k and 7 = e 1 ' 2 . (2) 
Define the height fluctuation field h^ nc (^,x) b}Q 



^We attempt to use capital letters for all variables (such as X, T) on the macroscopic level of the stochastic PDEs 
and polymers. Lower case letters (such as x, i) will denote WASEP variables, the microscopic discretization of these 
SPDEs. 
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Our first main result is a fluctuation theorem for the WASEP at the edge of the rarefaction fan. 

Theorem 3. Let p_ = 0, p+ = 1/2, and /i^ uc (^, x) be the edge fluctuation field defined in Definition 
Then for t,j and x as in (d)), and /i fluc as in ([3J) ; 

limP (h««%x) > 2-V3 ( x 2 - s j) = F^(s) 

where the edge crossover distribution F^^(s) is given by 

F°% e (s) = ^ e-» d j- det(I - ^ edge ) i2( r,). (4) 

The operator Kg d&e , which depends on T and X , and the contour f v , C are defined in Definition 
\18l Alternative formulas for the distribution function F^ d ^(s) are given in Section \5. S\ 



This theorem is proved in Section [3j The proof uses the same method as the proof of the main 
theorem of [2] and relies upon a recently discovered exact formula for the probability distribution 
for the location of a fixed particle in ASEP with step Bernoulli initial conditions (in our case p_ = 
and p+ = 1/2). The main technical modification is due to a new infinite product, which we call 
g(C)- Relating this to q— Gamma functions we are able to extract the new asymptotic kernel which 
now also contains Gamma functions. 

As we will see in Section 11.21 the F^y > (s) distribution is also the one-point distribution for the 
KPZ equation ([5]) with specific initial data (|15p . It is clear from this result that time and space 
scale differently. Specifically, the ratio of the scaling exponents for time : space : fluctuations is 
3:2:1. This scaling ratio was shown in [17J to hold for a wide class of 1 + 1 dimensional growth 
models. Finally, the X 2 shift with respect to the s variable reflects the parabolic curvature of the 
rarefaction fan nearby the edge. 

The above result should be compared (see Section 12. ip with the existing fluctuation theory for 
TASEP and ASEP. In those cases, using the same centering and scaling as in ([3]), [5J [9] and |42j 
(respectively for TASEP and ASEP) obtained formulas for the one-point probability distribution 
function. These formulas actually first arose in the study of the largest eigenvalue of rank one 
perturbations of complex Wishart random matrix ensembles [5]. Remarkably, the limiting distri- 
butions are the same regardless of the asymmetry 7, as long as it is held positive as the other 
variables scale to infinity. By scaling 7 as above we focus in on the crossover between the ASEP 
and the SSEP and the new family of edge crossover distribution functions represent this transition. 



For TASEP, [16J gave a formula for the asymptotic equal time height function fluctuation process 
(in terms of finite dimensional distributions). We paraphrase this as Theorem 1151 For our present 
case, it says that if we fix 7 = 1, m > 1 (the case m = 1 is just the [SJ [9] result mentioned above), 
then for any choices of T > 0, Xi, . . . , X m S R and s±, . . . , s m G R, 



m 



limP f){h^ c (±,x k ) > 2~ 1 I\XI - s k )} =P f|{^BM(X fc ) < s k } 



\fc=i / \fc=i 



where .A2-3.BM is a spatial process (defined below in Definition [20]) which interpolates between the 
Airy 2 process and Brownian motion. 
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Since in WASEP we scale the asymmetry with time and space in a critical way, the fluctuation 
distributions are not the same as for TASEP or ASEP. Rather than having 7 = e 1//2 one could 
perform asymptotics with 7 = ae 1 / 2 . Doing this is becomes apparent that increasing T is like in- 
creasing a, hence, one expects to recover the TASEP distributions from the WASEP edge crossover 
distributions as T /* 00. 

Conjecture 4. Let p_ = 0, p+ = 1/2, as well as t,^/,x and h^ uc (^,x) be as in Definition^ Then 
for any m > 1, 



T limJimP \T\{h^ c ^,x k ) > 2- 1 l\X 1 k - s k )} j = P f f| {A 2 ^Bu(X k ) < s k} 

where the joint distribution for the process ^2->bm is given in Definition\20\ in terms of a Fredholm 
determinant. 



Extracting asymptotics from the result of Theorem [3] we are able to confirm this conjecture in the 
case of m = 1 (see Section [5.11 for the proof): 

Corollary 5. The i*y v- e distribution has a long time limit which is given by 

hm F° d x e (s) = P (A 2 ->bm(X) < s) , 

where the above one-point function for the ^2^bm process coincides with the so called BBP- 
transition [5] in the study of perturbed Wishart random matrices and is given in Definition\2(A 

Thus by rescaling the edge crossover distribution for WASEP we recover the universal distribution 
at the edge of the rarefaction fan for TASEP and ASEP. In the other direction we can also extract 
the small T asymptotics, which are Gaussian. This is best stated in terms of the stochastic heat 
equation so we delay it to Proposition QTJ 



1.2. KPZ equation as the limit for WASEP height function fluctuations at the edge. 

Following the approach of [12], we prove that as e goes to zero, a slight variant on the fluctuation 
field h^ uc (-t,x) converges to the KPZ equation with appropriate initial data. 

The KPZ equation was introduced by Kardar-Parisi-Zhang in 1986 as arguably the simplest sto- 
chastic PDE which contained terms to account for the desired behavior of one-dimensional interface 
growth [25] : 

d T u = -\{d x nf + \d x n + w (5) 

where W(T,X) is space-time white noise (see [2] Section 1.4 for a rigorous definition of white noise) 

E[W(T, X)W(S, Y)} = 8{T - S)S(Y - X). 

Despite its simplicity, the KPZ equation has resisted analysis for quite some time. The reason is 
that, even for nice initial data, the solution at a later time T > will look locally like a Brownian 
motion in X. Hence the nonlinear term is ill-defined. 

In order to make sense of this KPZ equation we follow |12] and simply define the Hopf-Cole solution 
to the KPZ equation as 

H(T,X) = - log Z(T,X) (6) 
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where Z(T, X) is the well-defined [35] solution of the stochastic heat equation, 

8 T Z = \d\Z - ZW. (7) 

Starting ([5]) with initial data H(0, X) means starting ([7]) with initial data Z(0, X) = exp{— H(0, X)}. 
However, one is best advised not to think in terms of H for the initial data since here we will deal 
with initial data for Z (such as Dirac-delta functions) which do not have a well-defined logarithm. 

The stochastic partial differential equation ([7]) is shorthand for its integral version, 

/oo r-T roo 

p(T,X-X )Z{0,X)dX - / piT-^X-X^ZiT^X^WidXtdn), (8) 
-oo JO J— oo 

where p{T,X) = {2-nT)- 1 / 2 exp{-A" 2 /2T} is the heat kernel. Iterating, one obtains the chaos 
expansion (convergent in J2f 2 of the white noise W) 



oo 



Z(T,X)=J2(-l) n In(T,X) (9) 



n=0 



where 



„ „ n n 

I n {T,X):= / T\p{T i+1 -T i ,X i+l -X i )Z(0,X o )dX o T[W{dT i dX i ), (10) 

JA' n (T) f = t f = l 

A' n (T) := {(T u . . . ,T n ) : = T < T t < • • • < T n < T n+l = T}, 

and X n+ \ = X. 

1.2.1. Microscopic Hopf-Cole transform. We can now show that the WASEP height fluctuation 
field converges to KPZ, in the sense that its Hopf-Cole transform converges to the solution of the 
stochastic heat equation. This idea was first implemented for equilibrium initial conditions in |12j 
and is facilitated by the fact that the Hopf-Cole transform of the fluctuations (with lower order 
changes to the scalings) actually satisfies a discrete space, continuous time stochastic heat equation 
itself [20]. Specifically let 

u e = p + q-2^qp = i e+ i e 2 + ( e 3) (H) 
X e = \ log(g/p) = e 1 / 2 + ie 3 / 2 + 0(e 5 / 2 ), (12) 
where we recall that with asymmetry 7 = e 1 / 2 we must have q = \ + ^e 1//2 and p = \ — |e 1//2 . 

Define the random functions Z e (T,X) by setting 

Z e (T,X)=exp{-A e / l7 (^,e- 1 X) + ^^p}. (13) 

Since p_ = 0, p + = 1/2, /i 7 (0, x) is |x| for x < and a simple symmetric random walk for x > 0. 
Using that fact and the Taylor approximation for X e ~ e 1 / 2 we find that for X < 0, A e /i 7 (0, e~ 1 X) 
is like e~ 1 l 2 X and for X > it is converging to a standard Brownian motion B(X). Thus negating 
and exponentiating we see that Z e (0,X) converges to initial data lx>o ex P{ — B(X)}. 

Definition 6. The solution of KPZ with half-Brownian initial data is defined as 

H edge (T, X) := - log Z edse (T, X) (14) 
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where Z edgc (T, X) is the unique solution of the stochastic heat equation 

Q TZ edge = I^cdgc _ ^edge^ Z C<isC {0,X) = 1 X > exp{- B{X)} . (15) 

The formal initial conditions for the (equally formal) KPZ equation would be % edgc (0, X) = B(X) 
for X > and % edgc (0, X) = -oo for X < 0. 



Now observe that via the Taylor expansions of v e and X e we have 



This suggests that 



log Z e (T,X) = ^^(r^T^-iy) + "L + (l). 



lim ^uc ( eZ^T )e -l x) = W*»(T,X) - I 



T l/3 

To state this precisely, observe that the random functions Z e (T, X) above have discontinuities both 
in space and in time. If desired, one can linearly interpolate in space so that they become a jump 
process taking values in the space of continuous functions. But it does not really make things 
easier. The key point is that the jumps are small, so we use instead the space D u ([0, oo); D U (R)) 
where D refers to right continuous paths with left limits and D U (M) indicates that in space these 
functions are equipped with the topology of uniform convergence on compact sets. Let & e denote 
the probability measure on D u ([0, oo); D U (R)) corresponding to the process Z e (T, X). 

Theorem 7. £P e , e 6 (0, 1/4) are a tight family of measures and the unique limit point is supported 
on C([0, oo); C(R)) and corresponds to the solution of ^7ty with initial conditions i!5\). 

The proof of this theorem is a variation on that of [12] and [2] and is given in Section [4l 

Our second main result is a corollary of Theorems [3] and [7] and provides an exact formula for the 
one-point distributions of KPZ with half-Brownian initial data. 

Corollary 8. For each fixed T > 0, X £ R and s € R, 

where F^ x e (s) is given in Definition \1& As T / oo the above converges to P (A2->bm(X) < s) 
(see Definition [2U\) . 

Translating Conjecture [5] into the KPZ language gives 
Conjecture 9. For m > 1, T > 0, X\, . . . ,X m £ R and si,...,s m £ 1 
T KmP (n | ^ dsc (T,2V3rV3x fc )-|; ^ 2 _ 1/3(x , _ sfc) |^ = p ^ {A ^ MXk) < ^ 



We prove m = 1 as the second statement of Corollary [H 
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One expects [26J that for large y, 

1 - F^(y) ~ ce~h 3/ \ F^(-y) ~ «T< (16) 

So far, we have not able to obtain ()16[) from the determinental formula ([4]) for -Fj^ 6 - In fact, using 
only the determinental formulas, it is an open problem to show that i^ d ^ e (y) — > as y — > — oo. 
We only know it is true because of Corollary [8] together with ([8]) , fjlQf) , (|14h , which show that 
H ge (T,X) is a non-degenerate random variable. The problem is that unlike the Airy kernel used 
to define the Tracy- Widom distributions, the eigenvalues of our crossover kernels are not in [0, 1]. 
We can obtain some asymptotics at the other end. 

Proposition 10. There exists 01,02,03 < 00 such that for T > 1, 

1 - F^ c (2 1 / 3 y - 2 1 / 3 T- 1 /3l og 2) < ^T 1 / 2 ( e ~<^ 3/2 + e -^v). (17) 

The proof is in Section 15.41 The constants and dependence on T given here are not expected to 
be optimal. We remark that using the same methods as in our proof, one can compute the upper 
tail for f 6 and, with better constants, the upper tail for the distribution F^\ (which on recalls 
does not depend on X). 

In terms of moments of 

^edgc( T)2 l/3 T 2/3 X ) _ T 

from the convergence in law of the one point distribution, and the general lower semicontinuity, 
one obtains lower bounds 



lim inf E 

T /-oo 



^edge (r)2 l/3 r 2/3 X) _Ty 



> C P {X) > 



where C P {X) = E[(A2->bm(X)) p ]. The corresponding upper bounds do not come as easily, though 
presumably they could be derived by an appropriate asymptotic analysis of the Tracy- Widom 
formulas for ASEP. 

In the other direction we can also extract the small T asymptotics. Solving the regular heat equation 



(x-Xq) 2 



d T Z = \d\Z with initial data (USD gives / °° e JL. e- Bi - x ^dX { 



o- 



Proposition 11. As T \ 0, 

Z edge (T, X) = / e~ BiXo) dX + T^Z^iT, T~ l l 2 X) + o(T^ 4 ) 

Jo v2 , kT 

where, given B(X),X > 0, Z imt (T, X) is a Gaussian process with mean zero and covariance 

roo roo 

Cov(Z init (T,X),Z init (T,Y))= / 9(X,Y,X ,X' )dX dX' 

Jo Jo 

where 

e -4,{ x +y- x o-x Y ri e 4(i- s ) Ts — 

*(*,y, X oA) = ^ j o Vs{l _ s) ds. 
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The above proposition (proved in Section I5.5|) indicates that the fluctuations scale and behave 
differently at the two extremes of T /* oo and T \ 0. One could have started with an asymmetry 
of 7 = be 1 / 2 . It turns out that the effect on the limiting statistic of modulating this b term is the 
same as modulating T. Thus, as T goes to infinity, it is effectively like (up to an interchange of 
limits) increasing the asymmetry away from the weakly asymmetric range to the realm of positive 
asymmetry. This explains the T 1//3 and T 2 / 3 scaling of fluctuations and space in the large T limit. 
On the other hand, taking T to zero is like moving towards symmetry, and this explains the T 1 / 4 
and T 1//2 scalings of fluctuations and space in this limit. These two classes are called the KPZ and 
EW universality classes. Thus we see that the KPZ equation is, in fact, the universal mechanism 
for crossing between these classes. 

These results along with similar results of |12|.[8]. [2] and those contained in Section [2] below, 
provide overwhelming evidence that the Hopf-Cole solution ([6]) to the KPZ equation is the correct 
solution for modeling growth processes. There do exist other interpretations of the KPZ equation, 
but they all suffer from the fact that they lead to answers which do not yield the desired scaling 
properties and limit distributions [15] . 



Remark 12. There is a Feynman-Kac formula for the solution of the stochastic heat equation 
Z(T,X)=E TjX ' 



Z(0,X) :exp: j-jf W(t,b(t))dt 



where we make use of the Wick ordered exponential and where Er,x is the standard Wiener measure 
on b(t) ending at position X at time T. This partition function is rigorously defined in terms of the 
chaos expansion ([9]) or alternatively, as a limit of mollified versions of the white noise [10] . Hence 
Z(T,X) has an interpretation as a partition function, and —7i(T,X) an interpretation as a free 
energy, of a continuum directed random polymer. These can be shown to be universal limits of 
discrete polymers with rescaled temperature pQ. All of the results in this article have alternative 
and immediate interpretations in terms of these polymer models. This polymer perspective is more 
along the line of the approach taken in [2] . 



1.3. Applications to KPZ in equilibrium. We now consider the Hopf-Cole solution of KPZ 
corresponding with growth models with equilibrium or stationary initial conditions. The initial 
data for KPZ is given [12] by a two-sided Brownian motion B(X), I G I, or in other words, the 
stochastic heat equation is given initial data 

Z cq (0,A) = exp{-B(X)}. 

As always, B(X) is assumed independent of the space-time white noise W . Strictly speaking, this is 
not an equilibrium solution for KPZ because of global height shifts, but it is a genuine equilibrium 
[12j for the stochastic Burgers equation 

d T u = -^d x u 2 + x -d\u + 

formally satisfied by its derivative U(T,X) = dxH CCi {T,X) where 

H cq (T,X) = -logZ cq {T,X). 

In [8] it was shown that the variance of % eq (T, 0) is of the correct order. In particular, there are 
constants < C\ < C2 < oo such that 

CiT 2/3 < Var(^ e9 (T,0)) < C 2 T 2/3 . 
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At this time we are not able to obtain the distribution of H eq (T, X) because the corresponding 
formulas of Tracy and Widom [32] are not in the form of Fredholm determinants. However we will 
obtain some large deviation estimates and moment bounds. The idea is to represent Z eq in terms 
of solutions with half Brownian initial data, 

Z eq = Z + + Z_ 

where Z + and Z_ solve the stochastic heat equation with the same white noise and initial data 

Z±(0,X) = l xeR± exp{-B(X)}. 

Note that the two initial data are independent, and we have as an additional tool the following 
correlation inequality which is novel to our knowledge. At a heuristic level it is clear that any 
two increasing functions of white noise should be positively correlated. Using the Feynman-Kac 
(continuum polymer) interpretation of the stochastic heat equation it is physically clear that the 
solution is increasing in the white noise. 

Proposition 1. (FKG inequality for KPZ) Let Z\, Z2 be two solutions of the stochastic heat 
equation ([?]) with the same white noise W, but independent random initial data -2a (0, X), i? 2 (0, X). 
We make the technical assumption that the solution to the stochastic heat equation with initial 
data JJi(0,X) and Z 2 (0,X) can be approximated, in the sense of process-level convergence, by the 
rescaled exponential height functions U3\) for WASEP. Let%i(T,X) = — logZi(0, X) denote the 
corresponding Hopf-Cole solutions of KPZ. Then for any Ti,T 2 > 0, Xi,X 2 £ M and si,s 2 £ K, 

P(Zi(Ti,Xi) < Sl and ^ 2 (T 2 ,X 2 ) < s 2 ) > P{Z x {T x ,Xi) < s x )P{Z 2 {T 2 , X 2 ) < s 2 ) (18) 

where P denotes the probability with respect to the white noise as well as the initial data. In 
particular, we have 

P(fti(Ti,Xi) > si andU 2 {T 2 ,X 2 ) > s 2 ) > P(H 1 (T 1 ,X 1 ) > Sl )P{U 2 {T 2 , X 2 ) > s 2 ). (19) 

This proof uses the FKG inequality at the level of a discrete system which converges to the stochastic 
heat equation. We choose to use the WASEP approximation for the stochastic heat equation 
explained in this paper, though it would also be possible to prove this result via a discrete polymer 
approximation. The WASEP approximation assumption is not very restrictive. The work of |12} [2] 
and this paper show that a wide range of initial data fall into this class, and one should be able to 
expand this even more. We also remark that stronger forms of the above FKG inequality may be 
formulated and similarly proved, though we do not pursue this further here. 

Proof of Proposition [7J By assumption we can approximate the relevant solutions to the stochastic 
heat equation in terms of the WASEP as Z\ >t and Z 2 ^. The graphical construction of ASEP can 
be thought of as a priori setting an environment of attempted left and right jumps. However, 
for our purposes we think of first throwing a Poisson point process of attempted jumps and then 
assigning the jumps a direction (left or right) independently with probability q and p. There is a 
natural monotonicity in this construction which says that changing a right jump to a left jump will 
only increase the associated height function. Taking the approximations for the initial data to be 
independent of each other, this implies that the events Ai )t = {Zi t6 < Sj} are increasing events if 
one thinks of the Poisson process of attempted jumps as giving a (random) lattice and the jump 
directions as being 1 (left) or —1 (right). This jump lattice is infinite, however, with probability one 
only a finite portion of it affects the value of the two Zj ]e (Tj, X{). Therefore, with probability one 
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the FKG inequality applies to this setting because of the product structure of jump assignments 
on the attempted jump lattice. Since the Ai )t are increasing events, they are positively correlated. 
Taking the limit as e — > gives the desired result (fl~8j) at the continuum level. Since — log is a 
decreasing function, (|19p follows immediately from (|18p . □ 



1 - F|, d gc (2 1 / 3 y - 2 1 / 3 r~ 1 / 3 log 2)) 2 < P (n e(i (T, 0) - % < -T^y) (20) 



Proposition 13. For all y G K and T > 0, 

and 



< 2(l- F^ e (2^ 3 y - 2 1 /3t~ 1 / 3 log 2; 



(i^ dge (-2 1 / 3 y - 2 1 / 3 T- 1 / 3 log2)) 2 < P (n cq (T,0) - g > T^y) (21) 

< 2F^ d ge (-2 1 / 3 y-2 1 / 3 T- 1 /3i og 2). 

One can derive similar expressions to those above for other values of X 7^ 0, though presently we 
do not state such results. 

Corollary 14. There exist 01,02,03 < 00 such that for T > 1, 

P (n cq (T, 0) - I < -T 1/3 y) < ciT 1 / 2 ^"" 2 ^ 2 + e- CsTl/3y ). (22) 
Furthermore, for each p > 0, i/iere exists C p > suc/i i/iai /or sufficiently large T , 

e {(n e(i {T, 0) - iff] > c p r p / 3 . (23) 

Proof of Proposition \13l If 

n c *(T,X) = -logZ c *(T,X), n±(T,X) = -\ogZ ± {T,X), 

then using the increasing nature of the logarithm and the fact that 

2 mm(Z + , Z_ ) < Z + + Z_ < 2 max(Z + , Z_ ) 

we have the following simple, yet significant inequality which expresses the equilibrium solution to 
KPZ in terms of two coupled half Brownian solutions: 

-log2+min(H+,H_) < n cq < -log2 + max(^+,^_). (24) 

Thus 

p(V*(T,0)-|J >T 1/3 ?/) < p(max(^ + (r,0),^_(T,0))-| >T 1 / 3 y + log2) (25) 

< 2P(H + (T,0)-|>T 1 / 3 y + log2 N ). (26) 



In ([26]) we used that by symmetry, H + (T, 0) and %-iT, 0) have the same distribution, and the 
upper bound of (|20p follows by Corollary [5J For the lower bound in (|20|) we have, by (|19|) . 

P(V»(T, 0)-ff >T l ^y) > p(mm(H+(T,0),H-(T,0))-% > T 1 /^ + log 2) (27) 

> P + (T, 0) - g > r 1 / 3 ^ + log 2) ] 2 . (28) 
(|2ip is obtained from the lower bound of (|24p in exactly the same way. □ 
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Proof of Corollary \14\ The large deviation bound (|22p follows from (|20p and Proposition [TUJ To 
prove ([23]) . suppose that G and F are probability distribution functions satisfying 



1 - G(x) > (1 - F(c 2 x + c 3 )) 2 and G(-s) > F 2 (-c 2 x + c 3 ) 
for all x £ R for some ci , c 2 , C3 . We have the bound that 

/roo roo 
x p dG(x) =p x p ~ l {l - G(x) + G(-x))dx > 2 / x(l - F(c 2 x + c 3 )) 2 + F 2 (-c 2 x + c 3 )dx. 
Jo Jo 

Hence the right side of ([23]) is bounded below by 

P So" ^ 1 { (1 ~ F ^^ 2l/3x ~ 2l/3r_1/3l °g 2 )) 2 + (FtT(~ 21/3x ~ 2 1/3 T- 1 / 3 log2)) 2 | dx. 

By Fatou's lemma, the limit inferior as T 00 is greater than the same integral with the distribu- 
tion function replaced by the distribution function of -4, 2 -s.bm(0). Since the latter is strictly 
positive, this gives (|23l) . □ 



1.4. Outline. In the introduction we have focused on the WASEP with p_ = 0, p+ = 1/2 two- 
sided Bernoulli initial conditions and velocity v = so as to be at the edge of the rarefaction fan. 
For those parameters we described the height function fluctuations for WASEP, the link to the 
solution of the KPZ equation with specific initial data, and then the fluctuation theory for that 
solution. In Section [2] we explain the situation for general values of p_ < p+ and v either inside 
the rarefaction fan or at the edge (see Remark [T7|) . In Section [2] we explain how the connection 
between WASEP and KPZ generalizes as well. Section T2.2I contains all of the important definitions 
of kernels, contours and Airy-like processes used in the paper. Section [3] contains a heuristic and 
then rigorous proof of Theorem [3l Section U] contains a proof that the WASEP converges to the 
KPZ equation (the claimed results of Section fl.2p . Finally, Section [5] contains a proof of the large 
T asymptotics of the KPZ equation, as well as other tail and short time asymptotics of the edge 
crossover distributions. 
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2. The fluctuation characterization for two-sided Bernoulli WASEP and KPZ 



In the introduction we focused on a particular choice of two-sided Bernoulli initial conditions where 
p_ = and p+ = 1/2. By looking at velocity v = (which corresponds to the edge of the rarefaction 
fan) we uncovered a new family of edge crossover distributions and showed that the fluctuation 
process near the edge converges to the Hopf-Cole solution to the KPZ equation with half-Brownian 
initial data. 



CROSSOVER DISTRIBUTIONS AT THE EDGE OF THE RAREFACTION FAN 



13 



In this section we consider what happens for other choices of /?_ < p+ and v. Theorem [TBI the 
main result of this section, shows that under the same sort of scaling as present for TASEP ( [16] or 
Theorem 1151 below) . the WASEP height function fluctuations converge to three different crossover 
distributions: 



(1) The fan crossover distributions for fluctuations around velocities v G (2p- — 1,2/9+ — 1); 

(2) The edge crossover distributions for fluctuations around velocities v = 2p± — 1; 

(3) The equilibrium crossover distributions for />_ = p + = p and fluctuations around the char- 
acteristic v = 2p — 1. 



The three cases above also correspond to the three different possible KPZ limits of the WASEP 
height function fluctuations. The stochastic heat equation initial data in the three cases are 



(1) Z^(0,X) = 5 X=0 

(2) Z cd ^(0,X) = l x > exp{-B(X)} 

(3) Z^(0,X)=eM-B(X)}. 



We similarly label the associated continuum height function H(T, X) = — log Z(T, X) as rl (T, X), 
U edge {T,X) and rl c<i (T, X). 



2.1. Height function fluctuations for two-sided Bernoulli initial conditions. Before con- 
sidering the WASEP height function fluctuations for general values of p_ < p+ and v it is worth 
reviewing the analogous theory developed in [231 E31 US El EE] for the TASEP. For the TASEP, 
[29j conjectured (based on existing results [7] for the related PNG model) a characterization of the 
fluctuations of the height function for TASEP started with two-sided Bernoulli initial conditions. 
They identified three different limiting one-point distribution functions which depending on the 
region's location, i.e. whether the velocity v is chosen so as to be: (1) within the rarefaction fan; 
(2) at its edge; (3) in equilibrium (p_ = p+ = p), at the characteristic speed v = 2/9—1. This 
conjecture was proved in [191 [9]. An analogous theory for the multiple-point limit distribution 
structure was proved in [HI [16] and again only depended on region (l)-(3). 



Let us now paraphrase Theorem 2.1 of [16] (which also includes the main result of [6]). Define a 
positive velocity version of h^ uc (^,x) as 

hJ±,vt + x(l - v 2 ) 1 / 3 ) - ( ±&-t + x(l - v 2 ) 1 ^) 

j,flucft \ 2 L 

"-y,v V 7 > ■ ^ _ ^2)2/3^1/3 

The (1 — v 2 ) scaling term is the only new element and simply reflects the change of reference frame 
due to the non-zero velocity. When v = we recover h^ nc from (|3"j). 

Theorem 15 ([16] Theorem 2.1 (a), paraphrased). Fix 7 = 1 (TASEP), m > and p-,p+ € [0, 1] 
such that /9_ < p + . Then for any choices of T > 0, X\, . . . , X m £ R and s%,...,s m G R if we set 
t = e" 3 / 2 T and x k = 2 l /H 2 ^X k , 
(1) For v G (2p_ - 1,2/0+ - 1), 

limP (f){h^(t,x) > 2- 1 '\X 2 - s k )}) = P ((\{A 2 {X h ) < s k }) . (29) 



\k=l / \k=l 
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(2) For p_ < p + and v = 2p± — 1, 

limP (j\{h^,x) > 2- 1 /3 (X 2 _ Sfc)} ^J = p (j} {A ^ BM{ ± Xk) < Sk} ^j . (30) 

(3) For p^ = p + = p and v = 2p — 1, 
/ m \ / m \ 

limP ( f|{A*£(±,x) > 2_1/3 (^ 2 - ) = P ( fl^™ ^ + J • ( 31 ) 



vfc=l / \fc=l 



The definitions of the three limit processes are given in Definition [20] of Section ES The X\ in the 
^4 e q probability reflects a shift of that amount which was already present in the definition of that 
process originally given in [6] (in that paper this process was called -4 s tat) though for our purposes 
we find it more informative to call it ^4 e q)- 

For ASEP, where 7 is less than one but still strictly positive, significantly less is known rigorously. 
The only case where results analogous to Theorem 1151 have been proved is p_ = 0, m = 1, though, 
given those results of [21142], it is certainly reasonable to conjecture that Theorem [T5l holds for all 
7 > 0. 

As we saw in the introduction, there exists a critical scaling for 7 going to zero at which there arises 
new limits for the height function fluctuations which correspond to the Hopf-Cole solution to the 
KPZ equation. In that WASEP scaling we then have the following theorem, very much analogous 
to Theorem 1151 

Theorem 16. Fix 7 = e 1//2 and < p_ < p+ < 1. Then for any choices of T > 0, X € M and 
seRifwesett = <r?>l 2 T and x = 2 1 / 3 t 2 / 3 X, 

(1) For v e (2p_ - 1,2/)+ - 1), 

limP x) > 2~ 1 / 3 (L e +X 2 - s)) = Fl™ v2)2T>x (s) 

(see Definition^ where L e = 2- 1 /3^-i/3 \ og ( e -V^/ c ) W uh c = (2p + - 1 - v)~ l + (v - 2p_ + l) -1 . 
As T /• 00 we recover the right hand side of l2§\) , P {A2{Xj t ) < s). 

(2) For v = 2p±- 1, 

HmP (h^(±,x) > 2^ 3 (X 2 _ g) ) = F ^l 2)2Txis) 
(see Definition \TS\) . As T /* 00 we recover the right hand side of {W\), P (A2->bm(±X) < s). 

(3) For p- = p + = p and v = 2p — 1, 

HmP ,x) > 2-V 3 (X 2 - s)) = Fll v2)2TX (s) 

for which presently no formula exists (see however Section \1.3\ for some bounds). As T /*■ 00 we 
should recover the right hand side of \31\) . 



The case p- = 0, p + = 1 and v = above was previously solved in [2] and [321 [331 EH- It should 
be noted that F^x does not, in fact, depend on X. The logarithmic correction L e in the case of 
the fan is unique to that case and does not have a parallel in the analogous TASEP result. 
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Remark 17. The above theorem can be proved in two ways. For p_ = one can use the formulas 
of [12] directly and extract asymptotics (as we do for the case p + = 1/2 in Section [3]). Aside from 
some small technical modifications to the proof, the only change is that for v ^ one must center 
the asymptotic analysis around the point 

_ n_ + v\ _ 2x /2 

? \l-vj T(l-v) 2 ' 

In order to arrive at the full (i.e., also p_ > 0) theorem above, one cannot appeal to exact formulas 
because those that exist |43j are not in a form for which it is known how to do asymptotics. 
Instead we will prove that the height function fluctuation process converges to the solution to 
the KPZ equation with initial data depending only on the region (fan, edge or equilibrium along 
the characteristic). In the case of the fan and the edge, we have exact formulas for the one- 
point distributions for these KPZ solutions, in the case of the fan, obtained from the special case 
p- = 0, p + = 1, and in the case of the edge, obtained from the special case p~ = 0,p + = 1/2. 
The equilibrium result also follows in the same way, despite not having a formula for F eci . We will 
give the general velocity version of Theorem [7] in a forthcoming paper. Thus, strictly speaking, 
in this paper we only prove the above theorem (and likewise Proposition [1]) for (a) p_ = 0, and 
general p + and v (using asymptotic analysis of formulas of [12]); (b) general p_ and p + , yet v = 
(combining Theorem [7J or the analogous theorems of [2] and [12] , along with the exact statistics 
results determined herein or in [2] for the KPZ equation itself). 

2.2. Kernel and contour definitions. Here we collect the definitions of the kernels and contours 
used in the statement of the main results of this paper. 



Definition 18. The edge crossover distribution is defined as 

ce' 



F~f(s) = f e-^ det(J - iff dge ) 
Jc ft- 

The contour C is given as 

C = W }7T/2<e<37r/2 U {x ± i} x>0 
The contours T v , are given as 

f„ = {^+ir:rG(-oo J -l)U(l,oo)}ufJ (32) 

f c = {-^+ir :r£ (-oo,-l) U(l,oo)}ufJ, (33) 

where f f is a dimple which goes to the right of XT -1 / 3 and joins with the rest of the contour, and 

where is the same contour just shifted to the right by distance C3 (see Figure [2]). The constant 
C3 is defined henceforth as 

The kernel i^| dge acts on the function space L 2 (f r; ) through its kernel: 

Kfz%fj,fj')= (34) 

/ exp{ (C 3 - r/ 3 ) + sT 1/3 (C - r/)}2 1/3 / £ dt ) , /g s ^ ttt:^—- 

Jr c 3 1 ' ^ Di J e t_~ r 2 i/ 3 ^ _ 2 i/3XT-V3 C _ fj 
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Figure 2. The contours Tq and T v extend vertically along the lines C3/2 and — C3/2 
except for a dimple to the right of XT -1 / 3 so as to avoid the poles of the Gamma 
function (marked as solid black circles). 



or, equivalently, evaluating the t integral: 

K?*>(fj,ff) = (35) 

Jt c Pi 3 K V ) + Sl K Vn s in(vr2V3(c_^)) T(2^f,' - 2^XT~^) C ~ - fj' 

T(z) is the standard Gamma function, defined for Re(z) > by F(z) = J °° s z ~ 1 e~ s ds and extended 
by analytic continuation to C — {0, — 1, — 2, . . .}. To follow previous works we continue to use the 
letter gamma for contours, but always with a subscript to differentiate them from the Gamma 
function. 

Definition 19. The fan crossover distribution is defined as 

Ft^x(s) = J_e-^det{I-K { r) L ^t v y 
The kernel acts on the function space L 2 (T V ) through its kernel: 

K^V, fj') = J_ exp{-|(C 3 - rf ) + sT^(( ~ n')}^ ( /_°° ^"t'T^ * 
or, evaluating the inner integral, equivalently: 

*?■«.*)- / ex P{ -| g - e) +^«- 1)} * 2 :'^'*^' A 

Jr ( 3 sm(7r2 i /^(^ - rf)) ( - n 

As far as the choice of contours and T v , one can use the same as above without the extra 
dimple (since there are no poles to avoid now). This distribution is closely related to the crossover 



-00 



e 1 - » I C - V 
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K A2 (X,x;X',y) = < 



distribution of [2], though one observes that the scalings are slightly different. As we now have a 
whole class of crossover distributions we find it useful to name them more descriptively. 

Definition 20. The Airy2 process is defined in terms of finite dimensional distributions as 

P ^f]{A 2 {X k ) < s k }j = det(/-x*-KU 3 X»)x3({x 1 ,... P x- m }xR)> 
where Xs(X k ,x) = lr x>Sk \, and K^ 2 is the extended Airy kernel: 

' /-oo 

/ dte {x '- x)t Ai(t + x)Ai(t + y), X>X', 
< Jo ,o 

- dte ( - x '~ x)t Ai(t + x)Ai(t + y), X < X' . 

, J —oo 

The Airy2 process was discovered in the PNG model [28J. It is a stationary process with one-point 
distribution given by the GUE Tracy- Widom distribution .Fgue [38] ■ An integral representation of 
K_a 2 can be found in Proposition 2.3 of [23]; another form is in Definition 21 of [13] in the M = 
case. 

We denote by ^2^bm the transition process from Airy2 to Brownian Motion. It is defined in terms 
of finite dimensional distributions as 

p (^(~){-A2^BM(X k ) < s k }\ = det(J - XsK ABM ^ 2 Xs)L2({x 1 ,...,x m }xm.), 

where K A2 ^ BM is the rank-one perturbation K A2 : 

^bm (*> x; X',y) = K M (X, x; X',y) + Ai(x) (-e* x ' 3+x 'y - dtK\{t + y)e" x '^ . 

This transition process was derived in [22] and was shown to arise in TASEP at the edge of the 
rarefaction fan in [16J . An integral representation of the kernel can be found in [13], Definition 21, 
in the m = 1 case. For m = 1 this formula corresponds with the distribution F\ of [5]. 

The definition of the process for equilibrium TASEP, ^4 e qj is quite intricate and is given in [6] where 
it is actually called the .4 s tat process. Its joint distributions is the right hand side of equation (1.9) 
in [6]. 



3. Weakly asymmetric limit of the Tracy- Widom step Bernoulli ASEP formula 



In this section we will prove our main result, Theorem [3l The proof of that theorem follows by 
combining the proof of the main theorem of [2] (for step initial condition WASEP) with a few 
lemmas to cover a new element (the g{C) term stated below) which shows up for step Bernoulli 
initial conditions. As noted in the introduction, the key technical tool behind this proof is the exact 
formula for the transition probability of a single particle in ASEP with step Bernoulli initial data 
[42]. 

Theorem 21 (Main results of [32]). Let q > p with q+p = 1, 7 = q — p and r = p/q. Fix p_ = 
and for p + G (0, 1] set 

a = (1 - P+)/p+- 
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Since p_ = we can initially label our particles 1,2,3, .. . by setting the leftmost to be particle 1 
and the second left most to be particle 2, and so on. Let x(t, m) denote the location of particle m 
at time t. Then for m > 0, t > and i£Z, 1.2) gives the following exact formula 

P(x( 7 - 1 t,m)<x)= f ^n(l-/ir fc )det(J + M j;W,) (36) 

JS T+ V k = 

where S T + is a circle centered at zero of radius strictly between r and 1, and where the kernel of 
the determinant is given by 

4M)=[ exp{*(C) - *M}§^4^C (37) 
F Jr ( V [(-V) 9(0 

where 



' k 

z k 



k=—oo 

*(C) = A(C)-A(0 (38) 
A(C) = -xlog(l-C) + Y^ + ^logC 

g(o = n( i+rn «o- ( sg ) 

n=0 

The contours are a little tricky: ij and rf are on T v , a circle of diamete^ [—a~ l + 25, 1 — 5] for 5 
small. And the ( integral is onT^, a circle of diameter [— a^ 1 + q, 1 + <r] . One should choose q so as 
to ensure that \C/f]\ £ (l) r_1 )- This choice of contour avoids the poles of the new infinite product 
which are at —a~ l T~ n for n > 0. Of course we can take 5 to depend on e. 



3.1. Heuristic explanation of the asymptotics of the Tracy- Widom formula. We start 
by restating the result and give a heuristic explanation of the proof. In Section 13.21 we will give a 
complete proof of these asymptotics, roughly following the method of proof of Theorem 1 of [2]. 



The following theorem uses different scalings so as to conform to the notation of [2]. Therefore, the 
resulting formula differs and we introduce the kernel Ka SC ' ■ From the following result, by careful 
scaling, one arrives at Theorem [3l 

Theorem 22 (Equivalent to Theorem [3] after rescaling). Consider e > 0, T > and X G K and 

set p_ = 0, p+ = 1/2, t = e~ 3/2 T ; x = e^X and 7 = e 1 / 2 . Then 

djl 



limP (^[hy&x) - §] > -a) = limP(x 7 (i m) < x) 



c 



where a = a(s) = s + 



m 



-1/2 



X" 

-a H — — I H hx 

2T 2 



2 Following [35], this means the circle is symmetric about the real axis and intersects it at — a 1 +28 and 1 — 5. 
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and the operator Ka SC ' acts on the function space L 2 (f „) through its kernel: 

(40) 

The contours C, and T v are defined in Definition \18\ though for the last two contours, the 
dimples are modified to go to the right of the poles of the Gamma function above (the rightmost of 
which lies at l^^XjT). 

We now proceed with the heuristic proof of the above result. Note that given the values of p- and 
p+, the parameter a defined above in Theorem 1211 is equal to 1. We will, however, keep a in the 
calculations since one can then see readily how to generalize to a / 1. 

The first term in the integrand of (|36p is the infinite product J|^L (1 — pr k ). Observe that r = 
q/p ~ 1 — 2e 1 / 2 and that S T +, the contour on which fi lies, is a circle centered at zero of radius 
between r and 1. The infinite product is not well behaved along most of this contour, however we 
can deform the contour to one along which the product is not highly oscillatory. Care must be 
taken, however, since the Fredholm determinant has poles at every \x = r k . The deformation must 
avoid passing through them. As in [2] observe that if 

// = e 1/2 fi (41) 

then 

oo 

pr k ) « e -£~=o/- fe = e -M/(l-T) „ e -A/2. 

k=0 

We make the \x i— > e~ 1 / 2 jl change of variables and find that if we consider a jl contour 

C e = {e ?e } 7r / 2 <e<3 7 r/2 U {x ± i} 0<a; < e -l/a_l U {^ l/2 ~ 1 + %}-l <2/ <l, 

then the above approximations are reasonable. Thus the infinite product goes to exp{— fl/2}. 

Now we turn to the Fredholm determinant and determine a candidate for the pointwise limit of 
the kernel. The kernel Jj^(rj, if) is given by an integral whose integrand has four main components: 
An exponential 

exp{A(C)-A(r/)}, 

a rational function (we include the differential with this term for scaling purposes) 

d(/v'(( ~ V), 

a doubly infinite sum 
and an infinite product 

g(v')/g(C). 



We proceed by the method of steepest descent, so in order to determine the region along the ( 
and ij contours which affects the asymptotics we must consider the exponential term first. The 



20 



I. CORWIN AND J. QUASTEL 



argument of the exponential is given by A(£) — A(?/) where 

A(C) = -x log(l - C) + + m log(C), 

where x, t and m are as in Theoreml22i For small e, A(£) has a critical point in an e 1 / 2 neighborhood 
of -1. For purposes of having a nice ultimate answer we choose to center in on the point 

e = _l_ 2e V2£ 

We can rewrite the argument of the exponential as (A(£) — A(£)) — (A(?/) — A(£)) = ^Z>{() — ^(v')- 
The idea of extracting asymptotics for this term (which starts like those done in [H] but quickly 
becomes more involved due to the fact that r tends to 1 as e goes to zero) is then to deform the 
( and rj contours to lie along curves such that outside the scale e 1 / 2 around £, Re^C) is very 
negative, and Re^r/) is very positive, and hence the contribution from those parts of the contours 
is negligible. Rescaling around £ to blow up this e 1 / 2 scale, gives us the asymptotic exponential 
term. This change of variables sets the scale at which we should analyze the other three terms in 
the integrand for the J kernel. 



Returning to ^(C)) we make a Taylor expansion around £ and find that in a neighborhood of £ 

*(0--§ e - 3/2 (c-e) 3 + f e - 1/2 (c-c). 

This suggests the following change of variables 

C ~ = 2 -4/3 e -l/2 (c _ c) ~ = 2 -4/3 e -l/2 (r? _ ^ ~, = 2 -4/3 e -l/2 (r? / _ ^ (42) 

after which our Taylor expansion takes the form 

f(C)--|c 3 + 2 1 / 3 aC 

In the spirit of steepest descent analysis we would like the £ contour to leave £ in a direction where 
this Taylor expansion is decreasing rapidly. This is accomplished by leaving at an angle ±2tt/3. 
Likewise, since ^(r/) should increase rapidly, rj should leave £; at angle ±7r/3. Since p + = 1/2, a = 1 
which means that the £ contour is originally on a circle of diameter [—1 + 5, 1 + 5] and the r\ contour 
on a circle of diameter [— 1 + 25, 1 — 5] for some positive 5 (which can and should depend on e so as to 
ensure that \C/f]\ £ (l) r_1 ))- I n order to deform these contours to their steepest descent contours 
without changing the value of the determinant, great care must be taken to avoid the poles of /, 
which occur whenever (,/rf = r fe , k G Z, and the poles of 1/g, which occur whenever £ = — r~ n , 
n > 0. We will ignore these considerations in the formal calculation but will take them up more 
carefully in Section 13.21 The one very important consideration in this deformation, even formally, 
is that we must end up with contours which lie to the right of the poles of the 1/g function. 



Let us now assume that we can deform our contours to curves along which \& rapidly decays in £ 
and increases in rj, as we move along them away from £. If we apply the change of variables in (|42p 
the straight part of our contours become infinite rays at angles ±27r/3 and ±7r/3 which we call 
and r^. Note that this is not the actual definition of these contours which we use in the statement 
and proof of Theorem [3] because of the singularity problem mentioned above. 
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Applying this change of variables to the kernel of the Fredholm determinant changes the I? space 
and hence we must multiply the kernel by the Jacobian term 2 4//3 e 1 / 2 . We will include this term 
with the nf(fj,, z) term and take the e — > limit of that product. 

Before we consider that term, however, it is worth looking at the new infinite product term 
g(Tj')/g(C')- In order to do that let us consider the following. Set 

log a(c — xr) loga(c-yr) 

q=l-r, a= , b = . 

log q log q 

Then observe that 



TT l + q-QM-c + ar) = (g a ;g)oo = _ b _ a 

11 1 + a- r )n a (_ c + yr ) (qb-a)^ T a (a) [ Q> 



n=0 



= r g( 6 ) c (fe-a) losr = ri_ r (-r 1 log(ac) + c 1 y + o(r)) c ^i (y _ x)loRr+o(rloKr) ^ 
T q (a) ri_ r (-r^ 1 log(oc) + c~ l x + o(r)) 

where the q-Gamma function and the q-Pochhammer symbols are given by 
when \q\ < 1 and 

(a; q)oo = (1 - a)(l - aq)(l - aq 2 ) 

The notation o(/(r)) above refers to a function f'(r) such that f'(r)/f(r) — > as r — > 0. The 
q-Gamma function converges to the usual Gamma function as q — > 1, uniformly on compact sets 
(see [3] for more details and a statement of this result). 

Now consider the g terms and observe that in the rescaled variables this corresponds with (|43|) with 
r = 2e 1//2 , c = 1 (recall a = 1 as well) and 

Since ac = 1 and since we are away from the poles and zeros of the Gamma functions we find that 



iT" r(S"4) exp { 2l/3(e - f '' )log(2el/2 ^ (44> 

This exponential can be rewritten as 

exp | Z - log e | exp fe 1 / 3 log(2)(C - f?')} • (45) 

where 

~Z = 2^ 3 (C-V)- (46) 
It appears that there is a problem in these asymptotics as e goes to zero, however we will find 
that this apparent divergence exactly cancels with a similar term in the doubly infinite summation 
term asymptotics. We will now show how that loge in the exponent can be absorbed into the 
2^ e 1 / 2 nf(nXh') term. Recall 



k= 



k, 
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If we let no = Ll°g( e 1 ^ 2 )/ l°s( r )J then observe that 



k 

nf{ji,z)= y - ^' ,~ z k+n ° = z n °T n °fi y - — - h — 

fc=— oo k=— oo 



By the choice of no, T n ° ps e l l 2 so 

fif(n,z) ps z no flf(fi,z). 

The discussion on the exponential term indicates that it suffices to understand the behavior of this 
function only in the region where £ and rj are within a neighborhood of £ of order e 1 / 2 . Equivalently, 
letting z = C/V; it suffices to understand fj,f(n,z) ps z n °flf(fi,z) for 

_ c _ g + g/V^C ^ 1/2 - 

Let us now consider z n ° using the fact that logr ps — 2e 1 / 2 : 

z n ° ps (1 — e 1 / 2 ^) 6 " 172 ! 1 ^ 6 « e ~3 51 °s e . 

Plugging back in the value of z in terms of £ and j)' we see that this prefactor of z n ° exactly cancels 
the loge term which came from the g infinite product term. 

What remains is to determine the limit of 2 4/ ' 3 e 1/ ' 2 p>f{p>, z) as e goes to zero and for z ps 1 — e 1//2 z. 
This limit can be found by interpreting the infinite sum as a Riemann sum approximation for an 
appropriate integral. Define t = Zee 1 / 2 , then observe that 

00 uT te ~ 1/2 7 tt ~ 1 ' 2 , r°° rip- 2t p-~ zt 



-2t 



This used the fact that r te 1/2 — > e 2t and that z te 1/2 — >• e which hold at least pointwise in i. 
If we change variables of t to t/2 and multiply the top and bottom by e~* then we find that 

/oo ~ — zt/2 
-oo e - V 

As far as the final term, the rational expression, under the change of variables and zooming in on 
£, the factor of \jrj goes to -1 and the goes to j^=j- 

Therefore we formally find the following kernel: —K^ c, (fj,ff) acting on L 2 (f ?? ) where: 



r t I r°° r lf >-2 1/3 t(t-fi') \ r(2 1 / 3 C-f ) a? 



3 VS ,; vs \J-oo e*-/i j r(2V3^'-f)(_^' 

where a' = a + log 2 (recall that this log 2 came from (|45|) ) . 

We have the identity 



oo ~ -zt/2 

—. -dt = (-/i)- 2/2 7rcsc(7rz/2), (47) 

-oo e«-u 
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where the branch cut in Jl is along the positive real axis, hence (— /i)~ 2 / 2 = e — log ( — ^*' 2 where log 
is taken with the standard branch cut along the negative real axis. We may use the identity to 
rewrite the kernel as 



ft) = I exp{-|(C 3 - rf ) + 2 1 /3 a ' (C ~ - fj')}2^ 



(tt2 1 /3(^ - ff)) r (2V3^ _ *) g _ ~ ■ 



sm 



7T(-/i)- 2l/3 ^') r ( 2l/3 C-^j d(_ 

(48) 

To make this cleaner we replace /i/2 with /2. Taking into account this change of variables (it also 
changes the exp{— ji/2} in front of the determinant to exp{— fj,}), we find that our final answer is 

^e-^det(/-iC C '% ( f^ 
which, up to the definitions of the contours and Tq is the desired limiting formula. 

It is important to note the many possible pitfalls of such a heuristic computation: 1. Pointwise 
convergence of both the prefactor infinite product and the Fredholm determinant is not enough to 
prove convergence of the fi integral; 2. The deformations of the r) and £ contours to the steepest 
descent curves are invalid, as they pass through multiple poles of the kernel, coming both from 
the / term and the g term; 3. One has to show that the kernels converge in the sense of trace 
norm as opposed to just pointwise. The Riemann sum approximation argument can in fact be 
made rigorous though, in [2] an alternative proof of the validity of that limit is given via analysis 
of singularities and residues. 



These possible pitfalls are addressed below in Section [37 



3.2. Proof of Theorem [3]. In this section we provide a complete proof of Theorem [22l from 
which one recovers Theorem [3] via scaling. The proof follows the same argument to that of [2] . As 
a convention, c (or capitalized, primed etc. versions) will represent a finite constant which can vary 
line to line, unless explicitly noted. 

In Theorem [22j we have reformulated the claim of Theorem [3] in terms of the weakly asymmetric 
simple exclusion process with half step Bernoulli initial data. Our proof, therefore, reduces to a 
rigorous asymptotic analysis of Tracy and Widom's formula (|36|) . That formula contains an integral 
over a \i contour of a product of a prefactor infinite product and a Fredholm determinant. The first 
step towards taking the limit of this as e goes to zero is to control the prefactor, n^=o(-^ — A* 1 "^)- 
Initially [i lies on a contour S T + which is centered at zero and of radius between r and 1. Along 
this contour the partial products (i.e., product up to N) form a highly oscillatory sequence and 
hence it is hard to control the convergence of the sequence. 

The first step in our proof is to deform the fi contour S T + to the long, skinny cigar shaped contour 

C € = {e 1 / 2 e i0 } 7r/2 < e < 37r/2 U {x ± ie 1 / 2 } 0<;c < 1 _ el/2 U {1 - e 1 / 2 + e 1 / 2 iy}_ 1<2/<1 , 

(see Figure [3]). We orient C t counter-clockwise. Notice that this new contour still includes all of 
the poles at \i = r k associated with the / function in the J kernel. 

In order to justify replacing S T + by C e we need the following: 
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Figure 3. The S e contour is deformed to the C e contour via Cauchy 's theorem and 
then a change of variables leads to C e , with its infinite extension C. 

Lemma 23. In equation [36\) we can replace the contour S e with C e as the contour of integration 
for n without affecting the value of the integral. 

We thank the referee for pointing out a mistake in the proof of this result in [2], and suggesting an 
alternative proof which we detail in Section 13.4.21 

Having made this deformation of the /i contour, we now observe that the natural scale for /x is on 
order e 1//2 . With this in mind we make the change of variables 

Remark 24. Throughout the proof of this theorem and its lemmas and propositions, we will use 
the tilde to denote variables which are e 1 / 2 rescaled versions of the original, untilded variables. 



The fi variable now lives on the contour 

Ce = {e ie }7r/2<e<37r/2 U {x ± i} 0<x < £ -^-l U {e~ 1/2 ~ 1 + W}-l<y<l, 

which grow and ultimately approach 

C = W 9 }7T/2<e<37r/2 U {x ± i} x>0 - 

In order to show convergence of the integral as e goes to zero, we must consider two things, the 
convergence of the integrand for ft in some compact region around the origin on C, and the controlled 
decay of the integrand on C e outside of that compact region. This second consideration will allow 
us to approximate the integral by a finite integral in fi, while the first consideration will tell us what 
the limit of that integral is. When all is said and done, we will paste back in the remaining part of 
the fi integral and have our answer. With this in mind we give the following bound which is taken 
word for word from Lemma 2.3 of [2] and whose proof (given therein) relies only on elementary 
inequalities for the logarithm. 

Lemma 25. Define two regions, depending on a fixed parameter r > 1, 

Rl = {A = |Al < ■ // ln J 
sm(7r/10) 

R 2 = {ft: Re(/x) G [ f , e' 1 / 2 ], and Im(/7) G [-2,2]}. 
tan(7r/10) 
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Ri is compact and R\ U R2 contains all of the contour C t . Furthermore define the function (the 
infinite product after the change of variables) 

00 

Pe (A)=n(i-e i/ v fc ). 

k=0 

Then uniformly in fx G R\, 

p e ( M )-> e -A/2 (49) 

Also, for all e < eo (some positive constant) there exists a constant c such that for all fx G R2 we 
have the following tail bound: 

|p e (£)| < | e — A/2 J J e — ce^/ 2 A 3 J . (50) 

(By the choice of R2, for all fx G R2, Re(/z 2 ) > 5 > /or some /ixec? <5. TTie constant c can be taken 
to be 1/8.) 



We now turn our attention to the Fredholm determinant term in the integrand. Just as we did 
for the prefactor infinite product in Lemma [25] we must establish uniform convergence of the 
determinant for fx in a fixed compact region around the origin, and a suitable tail estimate valid 
outside that compact region. The tail estimate must be such that for each finite e, we can combine 
the two tail estimates (from the prefactor and from the determinant) and show that their integral 
over the tail part of C t is small and goes to zero as we enlarge the original compact region. For this 
we have the following two propositions (the first is the most substantial and is proved in Section 
while the second is proved in Section 13.4.2)) . 



Proposition 26. Fix s E R ; T > and left. Then for any compact subset of C we have that 
for all 8 > there exists an eo > such that for all e < eo and all fx in the compact subset, 

det(J + e^fxJ?^)^ - det(J - < S ' 

Here a' = a + log 2 and K^f c r is defined in equation |^0[ ) and depends implicitly on fx. 

Proposition 27. There exist c, d > and eo > such that for all e < cq and all fx G C e , 

p e {fx) det(J + e X/ Vfi^W,,) < c'e^ 



This exponential decay bound on the integrand shows that that, by choosing a suitably large (fixed) 
compact region around zero along the contour C t , it is possible to make the fx integral outside of 
this region arbitrarily small, uniformly in e G (0, eo). This means that we may assume henceforth 
that fx lies in a compact subset of C. 

Now that we are on a fixed compact set of fx, the first part of Lemma [251 and Proposition l26lcombine 
to show that the integrand converges uniformly to 

e~^ 2 r 

— _det(/-^ c ' ) L 2 ( f I)) 

and hence the integral converges to the integral with this integrand. 

To finish the proof of the limit in Theorem 1221 it is necessary that for any 5 we can find a suitably 
small eo such that the difference between the two sides of the limit differ by less than 5 for all e < eo . 
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Technically we are in the position of a 5/3 argument. One portion of 5/3 goes to the cost of cutting 
off the fl contour outside of some compact set. Another 5/3 goes to the uniform convergence of 
the integrand. The final portion goes to repairing the fx contour. As 5 gets smaller, the cut for 
the fl contour must occur further out. Therefore the limiting integral will be over the limit of the 
fl contours, which we called C. The final 5/3 is spent on the following Proposition, whose proof is 
given in Section [3.4.21 

Proposition 28. There exists c, d > such that for all fl £ C with \fl\ > 1, 

e -a/2 



(T„) 



< \c'e~ c n 



Recall that the kernel Ka SC,r is a function of fl. The argument used to prove this proposition 
immediately shows that Ka SC ' T is a trace class operator on L 2 (f r? ). 

It is an immediate corollary of this exponential tail bound that for sufficiently large compact sets 
of fl, the cost to include the rest of the fl contour is less than 5/3. This, along with the change of 
variables in fl described at the end of Section 13.11 finishes the proof of Theorem [22] 



3.3. Proof of Proposition 1261 In this section we provide all of the steps necessary to prove 
Proposition [26] To ease understanding of the argument we relegate more technical points to lemmas 
whose proof we delay to Section 13.4.31 

During the proof of this proposition, it is important to keep in mind that we are assuming that fl 
lies in a fixed compact subset of C. Recall that We proceed via the following strategy 

to find the limit of the Fredholm determinant as e goes to zero. The first step is to deform the 
contours and to suitable curves along which there exists a small region outside of which the 
kernel of our operator is exponentially small. This justifies cutting the contours off outside of this 
small region. We may then rescale everything so this small region becomes order one in size. Then 
we show uniform convergence of the kernel to the limiting kernel on the compact subset. Finally 
we need to show that we can complete the finite contour on which this limiting object is defined to 
an infinite contour without significantly changing the value of the determinant. 

Recall from Theorem [21] that is defined to be a circle of diameter [—1 + 1 + while is a 
circle of diameter [— 1 + 2?, 1 — The condition imposed on q is that for ( and rf on the above 
contours, \(/r/\ < (1,t _1 ). We take ? = e 1//2 /2, and since t _1 « 1 + 2e 1//2 , it is clear that for 
this choice of \C/f]\ < (1> T_1 )- The choice of contours is also such that the poles of the infinite 
product 1/(7(C), which occur at — r~ n for n > 0, lie to the left of the contours. Also recall 

s T 

The function fifJ-X/v') which shows up in the definition of the kernel for J has poles as every point 
(/r]' = z = T k for k G Z. 

As long as we simultaneously deform the contour as we deform so as to keep £/r/ away 
from these poles, we may use Proposition [42] (Proposition 1 of |41|). to justify the fact that the 
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determinant does not change under this deformation. In this way we may deform our contours to 
the following modified contours 1^,1^: 

Definition 29. Let r„^ and T^ i be two families (indexed by / > 0) of simple closed contours in C 
defined as follows. Let 

Both and Tri will be symmetric across the real axis, so we need only define them on the top 
half, r^/ begins on the real axis — 1 + e 1//2 and follows a smooth, northwesterly pointing curve and 
joins the vertical line with real part £ + e 1//2 /2 (see Figure [2] for an illustration of such a curve). It 
then follows the straight vertical line for a distance le 1 / 2 and then joins the curve 

\ + e l ' 2 {K{e)+a)\e ie (52) 

parametrized by 9 from 7r — le 1 / 2 + 0(e) to 0, and where a = — 1/2 + 0(e 1 ^ 2 ). The small errors 
are necessary to make sure that the curves join up at the end of the vertical section of the curve. 
We extend this to a closed contour by reflection through the real axis and orient it clockwise. We 
denote the first two parts (the northwesterly pointing curve and vertical line), of the contour by 
r^ r * and the remaining, roughly circular part by r^[ c . This means that T v j = U T^[ c , and 
along this contour we can think of parametrizing rj by 9 £ [0, it]. 

We define Y^ i similarly, except that it starts out at — 1 + e 1//2 /2, joins the vertical line with real 
part £ — e 1 ' 2 /2 and finally joins the curve given by equation (|52h where the value of 9 ranges from 
9 = 7r — /e 1//2 + O(e) to 9 = and where a = 1/2 + 0(e 1 ^ 2 ). We similarly denote this contour by 
the union of TfJ 1 and Tfj c . 



By virtue of these definitions, it is clear that e 1//2 |C/V — T k \ stays bounded away from zero for all 
k, that |C/?/| is bounded in an closed set contained in (1, r _1 ) for all £ £ F^i and rj £ r„ j, and that 
£ V 2 (C _|_ i) i s bounded from zero. Therefore, for any I > we may, by deforming both the rj and 
( contours simultaneously, assume that our operator acts on L 2 (r r?> ;) and that its kernel is defined 
via an integral along T^j. It is critical that we now show that, due to our choice of contours, we 
are able to forget about everything except for the northwesterly pointing curve and vertical part 
of the contours. To formulate this we have the following: 

Definition 30. Let xf er * an d xf rc be projection operators acting on L 2 (T rh i) which project onto 
L 2 (r^ rt ) and L 2 (r^"' c ) respectively. Also define two operators J^ ert ' r anc l jc irc > r wn i c h act on 
L 2 (r^^) and have kernels identical to J r (see equation (137]) ) except the ( integral is over r^ e z rt and 
-pare respectively. Thus we have a family (indexed by I > 0) of decompositions of our operator J 
as follows: 

jF _ T vert,Y vert , T vert,T arc , jcirc,T vert , T circ,T aire 
J — J l Xl + 'J I Xl + J l Xl T Jj Xl ■ 



We now show that it suffices to just consider the first part of this decomposition ( J^ er *' r x] ;er *) for 
sufficiently large I. 

Proposition 31. Assume that p, is restricted to a bounded subset of the contour C. For all 5 > 

there exist eo > and lo > such that for all e < eo and all I > Iq, 

| det(J + M^)L2(r, ;i! ) - det(J + J^ er ' ,r ) L 2 (r ,j t) | < 5. 
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Proof. As was explained in the introduction, if we let 

n = Llog(e- 1 / 2 )/log(r)J (53) 
then it follows from the invariance of the doubly infinite sum for /(//, z) that 

Hf^,z) = z n "(flf(il,z) + 0(e 1 / 2 )). (54) 
Note that the 0(e 1 ^ 2 ) does not play a significant role in what follows so we drop it. 

Using the above argument and the following three lemmas (which are proved in Section I3.4.3P we 
will be able to complete the proof of Proposition [3TJ 

Lemma 32. For all c > 0, there exist < Iq < oo and eo > such that for all I > Iq, e < eo and 

„ r- retire /■ c retire 

V fc 1 r),l ' s. fe 1 £,Z J 

Re(*(r/)) > c|e - v\e~ l/ \ Re(*(0) < ~c|£ - Ck" 1/2 - 

Additionally, there exists aC> 0, 0</*o<oo and cq > suc/i f/iai /or all I > Iq, e < eo and 

„ c ~pcirc /■ c pcirc 

V ^ L r/,l > ^ fc 1 £,Z ' 

Re(*(r?)) > C7|C - r/k" 1 , Re(tf(C)) < -C|£ - Ck" 1 . 

Lemma 33. For all I > there exist eo > and c > suc/i i/zai for all e < eo, rf G T^/ and 

IC-« I 

Lemma 34. For all c > i/iere exists Iq > and eo > suc/i i/iai /or a// Z > Zo, e < eo and a/*/* 
n G T^[ c and C € T^[ c 

Re ^n log(C/r / ) + £ log J < c(|£ - CI + K - "D^ 1 - 

The no above accounts for the z n ° from equation (|54p . Comparing the exponential (of order e _1 ) 
decay of the second part of Lemma [32] with the upper bound of Lemma [341 we find that since 
the constant in Lemma [341 is arbitrary, the decay of ^(n) overwhelms the possible growth of of 
(C/v) n ° di 7 ])/ div)- Additionally taking into account the polynomial control of Lemma [33] and the 
remaining term n), we find that for any 5 > 0, we can find Iq large enough that ||«^ ert ' r xf irc ||l> 

\\Ji irC,r xr r % an d ll^ CirC ' r xr rc |li are au bounded by 5/3. Technically, in order to show this we 
can factor these various operators into a product of Hilbert-Schmidt operators and then use the 
decay explained above to prove that each of the Hilbert-Schmidt norms goes to zero (for a similar 
argument, see the bottom of page 27 of [H]). This completes the proof of Proposition [31] □ 



We now return to the proof of Proposition [26] We have successfully restricted ourselves to con- 
sidering j^ er *' r acting on L 2 (r^ rt ). Having focused on the region of asymptotically non-trivial 
behavior, we can now rescale and show that the kernel converges to its limit, uniformly on the 
compact contour. 

Definition 35. Recall C3 = 2~ 4 / 3 and let 

n = £ + c 3 ~ V/ 2 ??, 7/ = $ + eg V/V, C = £ + C3 1 e 1 / 2 C 
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Under these change of variables the contours T^ rt and r^ ri become 

f = {c 3 /2 + ir:re {-c 3 l, 1) U (1, c 3 l)} U fj, 

f a = {-c 3 /2 + ir:re (-c 3 l, 1) U (1, c 3 l)} U t% 

where f ^ is a dimple which goes to the right of XT -1 / 3 and joins with the rest of the contour, and 
where is the same contour just shifted to the right by distance c 3 (see Figure [2|). 

As / increases to infinity, these contours approach their infinite versions, 

l\ = {c 3 /2 + ir:re (-oo, -1) U (1, oo)} U f J, 
f c = {-c 3 /2 + ir:r£ (-oo, -1) U (1, oo)} U fj. 
With respect to the change of variables define an operator jf acting on L 2 (T V ) via the kernel: 

U f( u W l£l/2 ^ 

Jr cl (£ + c 3 e L / 2 r]'){C-v) 

Lastly, define the operator xi which projects L 2 (T V ) onto L 2 (r v j). 

It is clear that applying the change of variables, the Fredholm determinant det(/ + Ji ert ' T )L 2 {T v ^ t ) 



becomes det(J + Xiii3\ Xi)l 2 (t 



We now state a proposition which gives, with respect to these fixed contours f„ / and j, the limit 
of the determinant in terms of the uniform limit of the kernel. Since all contours in question are 
finite, uniform convergence of the kernel suffices to show trace class convergence of the operators 
and hence convergence of the determinant. 

Recall the definition of the operator K^ C ' T given in equation fl40f) . For the purposes of this propo- 
sition, modify the kernel so that the integration in £ occurs now only over T^j and not all of T^. 
Call this modified operator . 

Proposition 36. For all 5 > there exist eo > and Iq > such that for all e < eo, I > lo, and jl 
in our fixed compact subset ofC, 



det(J + xi^Ji Xl) L 2(r V:1 ) ~ det(J - XlKZ' X'W ,) 



where a' = a + log 2. 



Proof. The proof of this proposition relies on showing the uniform convergence of the kernel of u jf 
to the kernel of K C ^ T , which suffices because of the compact contour. Furthermore, since the Q 
integration is itself over a compact set, it suffices to show uniform convergence of this integrand. 
The two lemmas stated below will imply such uniform convergence and hence complete this proof. 



First, however, recall that fJ,f(fJ,, z) = z n °(flf(fl, z) + 0(e 1 ^ 2 )) where no is defined in equation ([53 
We are interested in having z = C/ 7 /) which, under the change of variables can be written as 

z = 1 - e i/a 5 + 0(e), z = \C ~ ff) = 2 4 / 3 (C - if). 
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Therefore, since no = — ilog(e 1//2 )e x l 2 + 0(1) it follows that 



exp{-2 1 /3(c - ff) logfe-^JX! + G (l)). 



This expansion still contains an e and hence the argument blows up as e goes to zero. This exactly 
counteracts the asymptotics of the ratio of g{r}')/g{C) due to the following, 

Lemma 37. For all I > and all 8 > there exists eo > such that for all ff E r^z and £ E j 
we /ioue /or < e < eo, 



exp{n MC/r, )}-^ - r(2l/3 ^_x^ 



This is combined with the following two lemmas, and all three are proved in Section [3.4,31 

Lemma 38. For all I > and all 8 > there exists eo > such that for all rj' E T^z and £ E T^ 
we /iat> e /or < e < eo, 



*(C) - 



( C3 _^3 ) + 2 l/, a(c _~ ) 



1/3, 



< J. 



Similarly we have 



e *(C)-*0O _ e -i(C 3 -^ 3 )+ 2l/3a (C-*7') 



< 5. 



Lemma 39. For all I > and all 8 > f/iere exzsis eo > suc/i i/iai /or aZZ 77' E f ^ and £ E f £ j 
we have for < e < eo, 



e 1/2 /i/ I A 



00 fle- 2l/3t ^-^ 
-00 A 



< 8. 



The integral above in Lemma [39] converges since our choices of contours for £ and 77' ensure that 
Re(— 2 1//3 (£ — r}')) = 1/2. Note that the above integral also has a representation ()47j) in terms of 
the cosecant function. This provides the analytic extension of the integral to all z £ 2Z where 
z = 2 4 / 3 (£ — rj') (note, however, that we do not require use of this analytic extension due to our 
choice of contours). 



Finally, the sign change in front of the kernel of the Fredholm determinant comes from the 1/rj 1 
term which, under the change of variables converges uniformly to —1. The reason why a' = a + log 2 
arises here is due to the term e 2l/3 lo s( 2 )(C~'? ) from Lemma [371 This proves the desired result. □ 



Having successfully taken the e to zero limit, all that now remains is to paste the rest of the 
contours and to their abbreviated versions T v j and T^. To justify this we must show that 
the inclusion of the rest of these contours does not significantly affect the Fredholm determinant. 
Just as in the proof of Proposition [31] we have three operators which we must re- include at provably 
small cost. Each of these operators, however, can be factored into the product of Hilbert Schmidt 
operators and then an analysis similar to that done following Lemma [33] (see in particular page 
27-28 of [H]) shows that because Re(£ 3 ) grows like |C| 2 along (and likewise but opposite for n') 
there is sufficiently strong exponential decay to show that the trace norms of these three additional 
kernels can be made arbitrarily small by taking I large enough. 
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This last estimate completes the proof of Proposition [26l 
3.4. Technical lemmas, propositions and proofs. 

3.4.1. Properties of Fredholm determinants. Before beginning the proofs of the propositions and 
lemmas, we give the definitions and some important properties for Fredholm determinants, trace 
class operators and Hilbert-Schmidt operators. For a more complete treatment of this theory see, 
for example, [36] . 

Consider a (separable) Hilbert space % with bounded linear operators C{T~L). If A G C(H), let 
\A\ = V A* A be the unique positive square-root. We say that A G Bi(H), the trace class operators, 
if the trace norm ||^4||i < oo. Recall that this norm is defined relative to an orthonormal basis of 
H as \\A\\i := X]^Li( e n,; |A|e n ). This norm is well defined and does not depend on the choice of 
orthonormal basis {e„} n >i. For A G £>i(%), one can then define the trace tic A := YH?=i( e ni Ae n ). 
We say that A G B2(rl), the Hilbert-Schmidt operators, if the Hilbert-Schmidt norm \\A\\2 := 
Vtr(|Jl| 2 ) < oo. 

For A G Bi(Ti) we can also define a Fredholm determinant det(I + A)u. Consider Ui £ H and 
define the tensor product U\ <g) • • • <g) u n by its action on v\ , . . . , v n G % as 

n 

Ul <g> • • • <g> U n (vi, ... 1 V n ) = 

i=l 

Then <S>7=i % * s * ne s P an of all such tensor products. There is a vector subspace of this space which 
is known as the alternating product: 

n n 

/\{-H) = {h G (g)^ : Va G S n ,ah = -h}, 

i=l 

where au\ (g) • • • <g) u n = u a n\ <g> • • • (g) ii CT ( n ). If e±, . . . , e n is a basis for H then A • • • A ei k for 
1 < ii < . . . < ife < n form a basis of /\ n (H). Given an operator A G define 

r n (A)(ni ® • • • <8> «n) := ^4«i (8) • • • <8> Au n . 

Note that any element in /\ n (%) can be written as an antisymmetrization of tensor products. 
Then it follows that T n (A) restricts to an operator from /\ n (H) into /\ n (H)- If A G B\(H), then 
trT^(A) < ll^HJ/n!, and we can define 

oo 

det(/+ A) = l + ^tr(r(*)(A)). 
fc=i 

As one expects, det(/+ A) = Ylj(l + Xj) where Aj are the eigenvalues of A counted with algebraic 
multiplicity (Thm XIII. 106, [30]). 

Lemma 40 (Ch. 3 |36J). A i— > det(/ + A) is a continuous function on B\(T-L). Explicitly, 

|det(J + A)-det(J + B)| < \\A - B\\ x exp(||A||i + \\B\\ X + 1). 
If A G Bi(H) and A = BC with B,C G B 2 (W) tfien 

1 1-4| 1 1 < ||£||2||C|| 2 . 
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For A £ B\(J-L), 

|det(J + A)| < e m \ 

If A £ B2CH) with kernel A(x,y) then 

■ ' 2 

\A(x,y)\' 2 dxdy 

Lemma 41. If K is an operator acting on a contour S and x is a projection operator unto a 
subinterval of S then 

det(I + ^x)l2(s, m ) = det(J + x^x)l2(s, m )- 



In performing steepest descent analysis on Fredholm determinants, the following proposition allows 
one to deform contours to descent curves. 

Lemma 42 (Proposition 1 of [H]). Suppose s — ?■ T s is a deformation of closed curves and a kernel 
L(j],r)') is analytic in a neighborhood ofT s x T s C C 2 for each s. Then the Fredholm determinant 
of L acting on T s is independent of s. 



3.4.2. Proofs from Section \3.2\ We now turn to the proofs of the previously stated lemmas and 
propsitions. 



Proof of Lemma\23\ We thank the referee for suggesting the following simple proof of this result. 
The lemma follows from Cauchy's theorem once we show that for fixed e, the expression n^=o(^ — 
fj,r k ) det(I + /U J^) is analytic in \i between S e and C € (note that we now include a subscript \i on J 
to emphasize the dependence of the kernel on However, this expression was derived from and 
is equivalent to 

det(/ - \K) 



where A = T~ m \x and K is the operator (1) of [32]. The operator K does not depend on A and 
is thus is an entire function of A. Therefore the only singularities are at A = 1, . . . ,r~ m+1 , which 
correspond to [i = r, . . . ,r m . None of these singularities are between the two contours, thus the 
desired result follows. □ 



Proof of Lemma \2h\ We prove this with the scaling parameter r = 1 as the general case follows in 
a similar way. Consider 

00 

log( 5e ( / i)) = ^log(l-e 1 /V e fc ). 

k=0 
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We have J2T=o e x l 2 T k = ±(1 + e l / 2 c e ) where c e = 0(1). So for /2 G #i we have 

oo 



log(^)) + |(l + e 1/2 c e 



fc=0 

oo 



< ^llogCl-eV^ + eVV 



< 



1/2 - ■ k\2 _ e lF-l 



fc=0 

oo 

fc=0 

ce 1 / ~ | /J- 1 ~ See 



1 - r 2 4 - 4eV2 



< ce^lul 2 < cV/ 2 . 



The second inequality uses the fact that for \z\ < 1/2, | log(l — z) + z| < |z| 2 . Since jl G i?i it 
follows that |z| = e 1//2 |/i| is bounded by 1/2 for small enough e. The constants here are finite and 
do not depend on any of the parameters. This proves equation (149D and shows that the convergence 
is uniform in jl on R\. 



We now turn to the second inequality, equation (|50p . Consider the region, 

D = {z : arg(z) G [-^, £]} n : Im(z) G (-i, i)} n : Re(z) < 1}. 

For all z £ D, 

Re(log(l - z)) < Re(-z - z 2 /2). 
For fl G i?2j it is clear that e 1//2 /i G D. Therefore, using (|55|) . 

oo 

Re(log( 5e (A))) = ^RellogCl-e 1 / 2 ^)] 

fc=0 

oo 

< Y, (-Reie 1 / 2 ^} - Re[(e 1 / 2 /ir fc ) 2 /2] 



(55) 



fc=0 



1 1/2 



< -Re^/2) - -e^Re(^). 

This proves equation (|50() . Note that from the definition of R2 we can calculate the argument of fl 
and we see that | arg/2| < arctan(2tan(jg)) < j and \ jl\ > r > 1. Therefore Re(/i 2 ) is positive and 
bounded away from zero for all jl G i?2- □ 



Proof of Proposition 21. This proof proceeds in a similar manner to the proof of Proposition! 
however, since in this case we have to deal with e going to zero and changing contours, it is, by 
necessity, a little more complicated. For this reason we encourage readers to first study the simpler 
proof of Proposition [28l 

In that proof we factor our operator into two pieces. Then, using the decay of the exponential 
term, and the control over the size of the esc term, we are able to show that the Hilbert-Schmidt 
norm of the first factor is finite and that for the second factor it is bounded by |/i| Q for a < 1 (we 
show it for a = 1/2 though any a > works, just with constant getting large as a \ 0). This gives 
an estimate on the trace norm of the operator, which, by exponentiating, gives an upper bound 
e c \fi\ a on the size of the determinant. This upper bound is beat by the exponential decay in jl of 
the prefactor term p e . 
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For the proof of Proposition 1271 we do the same sort of factorization of our operator into AB, where 
here, 



A(C,r,) 



C-v 

with < d < 1 fixed, and 

B(r,,0 = e ^*®e*M-*WtfMv)^. 

We must be careful in keeping track of the contours on which these operators act. It will be 
convenient for this proof to move the contours for rj and £ to contours V® t and t which are 
defined in the same manner as T^j and T^ i in Definition [29] The difference, however, is that these 
new contours start at —1 + e 1//2 (resp. — 1 + e 1 / 2 /2) and go straight up for distance Ze 1//2 before 
joining k{0) with a = —1/2 + 0(e 1 / 2 ) (resp. a = 1/2 + 0{e 1 ^ 2 )). The purpose of this is to avoid 
the necessity of creating a dimple in the contours, thus allowing us to apply Lemma H3l in the form 
it was stated in [2\. Changing the location of this vertical portion of our contours does not affect 
the Taylor expansion we performed since we can still center our rescaled variables at £, and all of 
those results were valid in a compact region with respect to the rescaled variables. 

Now using the estimates of Lemmas 1321 and 1381 we compute that ||^4||2 < oo (uniformly in e < eo 
and, trivially, also in jj). Here we calculate the Hilbert-Schmidt norm using Lemma EEUl Intuitively 
this norm is uniformly bounded as e goes to zero, because, while the denominator blows up as 
badly as e -1 / 2 , the numerator is roughly supported only on a region of measure e 1//2 (owing to the 
exponential decay of the exponential when £ differs from £ by more than order e 1 / 2 ). 

We wish to control ||i?||2 now. Using the discussion before Lemma [32] we may rewrite B as 

\vj 9{C) V 

Lemmas [32] and [38] show that 

| e -c*(O e *(0-*(ij)| < e-e-^df-ijI+K-CI) 
for some fixed constant C > 0. On the other hand, Lemma [34] shows that 

C\ no aW)_ ] < e _ e -i C '(|{-^[+ie-c[) 



vj g{() 

for any constant C . In particular, we can take C < C which shows that all of the term (besides 
the / term) decay exponentially fast. 

The final ingredient in proving our proposition is, therefore, control of \fif(fl,z)\ for z = C/v'- We 
break it up into two regions of r/, £: The first (1) when \r}' — C\ < c for a very small constant c and 
the second (2) when |r/ — Q\ > c. We will compute ||-B||2 as the square root of 

/ \B( V ,()\ 2 d7 1 d(+ f \B( V ,()\ 2 dr]d(. (56) 

Jr)X&Case (1) J»y,CeCase (2) 

We will show that the first term can be bounded by d'\jl\ 2c " for any a < 1, while the second term 
can be bounded by a large constant. As a result ||-E>||2 < d'\jl\ a which is exactly as desired since 
then ||AB||i < e c "l^ a ( see Lemma l40l) . 
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Consider case (1) where \rf — Q\ < c for a constant c which is positive but small (depending on 
T). One may easily check from the definition of the contours that e _1//2 (|C/ 7 ?l — 1) is contained 
in a compact subset of (0,2). In fact, Q/n' almost exactly lies along the curve \z\ = 1 + e 1 / 2 and 
in particular (by taking eo and c small enough) we can assume that Q/n never leaves the region 
bounded by \z\ = 1 + (1 ± r)e 1//2 for any fixed r < 1. Let us call this region R t ^ r . Then we have 

Lemma 43. Fix eo and r £ (0, 1). Then for all e < eo, fi € C e and z £ R e , r , 

<c\p\ a /\\-z\ 

for some a G (0, 1), with c = c(a) independent of z, fx and e. 

Remark 44. By changing the value of a in the definition of k{6) (which then goes into the definition 
of l and T° ; ) and also focusing the region R 6 ^ r around \z\ = 1 + 2ae 1 / 2 , we can take a arbitrarily 
small in the above lemma at a cost of increasing the constant c = c(a) (the same also applies for 

Proposition I28p . The \ fi\ a comes from the fact that (1 + 2ae 1 / 2 ) 2 e 1 log l^l « \ f\ a . Another remark 
is that the proof below can be used to provide an alternative proof of Lemma [39] by studying the 
convergence of the Riemann sum directly rather than by using functional equation properties of / 
and the analytic continuations. 



We complete the ongoing proof of Proposition [27] and then return to the proof of the above lemma. 

Case (1) is now done since we can estimate the first integral in equation ([55]) using Lemma 03] and 
the exponential decay of the exponential term outside of \rj — C\ = 0(e 1 ^ 2 ). Therefore, just as with 
the A operator, the e -1 / 2 blowup of \fif{Jl, C/v')\ is countered by the decay of the other terms and 
we are just left with a large constant time \ fi\ a . 

Turning to case (2) we need to show that the second integral in equation (I56p is bounded uniformly 
in e and l~l G C e . This case corresponds to \rj' — (\ > c for some fixed but small constant c. Since 
e 1//2 (|C/ 7 ?I — 1) stays bounded in a compact set, using an argument almost identical to the proof of 
Lemma[33]we can show that \fif{fi, C/v)\ can be bounded by C|/i| c ' for positive yet finite constants 
C and C . The important point here is that there is only a finite power of |/2|. Since \fi\ < e -1 / 2 
this means that this term can blow up at most polynomially in e" 1 / 2 . On the other hand we know 
that the combination of the other terms decay exponentially fast like e~ e c for some small yet 
finite constant c. Hence the second integral in equation (156p goes to zero. 

We now return to the proof of Lemma 03] which will complete the proof of Proposition [271 



Proof of Lemma \43[ We will prove the desired estimate for z with \z\ = 1 + e 1//2 . The proof for 
general z £ R t>r follows similarly. Recall that 



We will recenter doubly infinite sum at around the value 

fc* = LK 1 / 2 io g HJ. 
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This is motivated by the fact that for Im/2 = dbl, and large real part, the denominator is approxi- 
mately minimized when r fc = l/|/x|, corresponding to k ~ k*. This centering results in 



By the definition of k*, 
Thus we find that 

where 



\ Z f = | A |1/2 (1 + 0(e l/ 2)) 

|Af(MI = l/i| 1/2 EE 



-oo 1 — cjr 



and is roughly on the unit circle except for a small dimple near 1. To be more precise, due to the 
rounding in the definition of k* the oj is not exactly on the unit circle, however we do have 

>e 1 / 2 , \u\-l = 0{e 1 ' 2 ). 

The section of C e in which p, = e -1 / 2 — 1 + iy for y £ (—1, 1) corresponds to oj lying along a small 
dimple around 1 (and still respects |1 — oj\ > e 1 / 2 ). We call the curve on which u lies fl. 

We can bring the \p\ 1 ^ 2 factor to the left and split the summation into three parts, so that 



I Ml 



-1/2 



|£/0M)l equals 



-OO l-uJT k ^fc: 



'1/2 



fc , V^oo orr* 1 Jl 



We will control each of these term separately. The first and the third are easiest. Consider 



(*-i)Ei 



-1/2 



-OO l- WT fc ' 

We wish to show this is bounded by a constant which is independent of ft and e. Summing by parts 
the argument of the absolute value can be written as 

,„-r-e" 1/2 +l -,-1/2+1 , f-1 \V^-f" 1/2 t- 



1 — LOT 



-e-l/2+l' 



/fc = -00 (l- Wr fe )(l-WT fe + l)' 



We have r _<E 1/2+1 e 2 and \z~ e L,z+1 \ e _1 (where e ~ 2.718). The denominator of the first 
term is therefore bounded away from zero. Thus the absolute value of this term is bounded by a 
constant. For the second term of (157p we can bring the absolute value inside the summation to get 



|l-r|££- 



1/2 



(l-LOT k )(l-UJT k + 1 ) 



The term 



LOT 

{l-LOT k ){l-LOT k + 1 ) 



stays bounded above by a constant times the value at k 



-1/2 



Therefore, replacing this by a constant, we can sum in \z\ and we get ^~ito ■ The numerator, 

as noted before, is like e _1 but the denominator is like e 1//2 /2. This is cancelled by the term 
1 — r = 0(e 1 ^ 2 ) in front. Thus the absolute value is bounded. 



- e -V2 



The argument for the third term of equation (I57D works in the same way, except rather than 
multiplying by |1 — z\ and showing the result is constant, we multiply by |1 — rz\. This is, however, 
sufficient since |1 — rz\ and |1 — z\ are effectively the same for z near 1 which is where our desired 
bound must be shown carefully. 
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We now turn to the middle term in equation (J57j), which is more difficult. We will show that 



:i-^)EL'-./^^ =o(iog|AD- 



This is of smaller order than \p,\ raised to any positive real power and thus finishes the proof. For 
the sake of simplicity we will first show this with z = 1 + e 1 / 2 . The general argument for points z 
of the same radius and non-zero angle is very similar as we will observe at the end of the proof. 
For the special choice of z, the prefactor 1 — z = e 1 / 2 . 

The idea, as mentioned in the heuristic proof, is to show that this sum is well approximated by a 
Riemann sum. In fact, the argument below can be used to make that formal observation rigorous, 
and thus provides an alternative method to the complex analytic approach we take in the proof of 
Lemma [321 The sum we have is given by 

f l/2V £ - 1/2 cjt" fc _ 1/2 Y^-V 2 w{\-eW+Q{e))* ( , ?) 

6 2^fe=-e-V2 ~ e 2^fc=-e-V2 i_ w (l-2eV2+0( £ )) fc ^ ' > 

where we have used the fact that tz = 1 — e 1 / 2 + 0(e). If k = te~ 1 / 2 then the sum is close to a 
Riemann sum for 

S\ T^dt. (58) 

We use this formal relationship to prove that the sum in equation (|57|) is 0(log|/2|). We do this 
in a few steps. The first step is to consider the difference between each term in our sum and the 
analogous term in a Riemann sum for the integral. After estimating the difference we show that 
this can be summed over k and gives us a finite error. The second step is to estimate the error of 
this Riemann sum approximation to the actual integral. The final step is to note that 

£ lT ^F^~|log(l-a;)|~log|A| 

for ui G f2 (in particular where |1 — ui\ > e 1 / 2 ). Hence it is easy to check that it is smaller than any 
power of |/2|. 

A single term in the Riemann sum for the integral looks like e 1 / 2 ue~ ktl/2 / '(1 — ue~ 2kel/2 ) . Thus 
we are interested in estimating 



,1/2 



w(l-e 1 / 2 +Q(e)) k _ u>e- k ^' 2 
l-w(l-2e 1 /2+0(e))fe l-uje- 2k * 1/2 



(59) 



We claim that there exists C < oo, independent of e and k satisfying fee 1 / 2 < 1, such that the 
previous line is bounded above by 



(l-cj+a>2fceV2) (l-ui+^ke 1 / 2 ) 2 ' 

To prove that (|59l) < (!60p we expand the powers of k and the exponentials. For the numerator and 
denominator of the first term inside of the absolute value in (|59p we have w(l - e 1 / 2 + 0{e)) k = 
uj - ujke 1/2 + 0(k 2 e) and 



l-w(l-2e 1/2 +0(e)) fe = 1 - u + oj2ke 1/2 - u2k 2 e + 0{ke) + 0(k 3 e 3/2 ) 

e+Q(kt)+0(k 3 
l-uj+^ke 1 / 2 



[l-u; + u2ke l ' 2 )(l - " 2fc2£+ °^ + °( fc3£3/2 >- 
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Using 1/(1 — z) = 1 + z + 0(z ) for |z| < 1 we see that 

cj(l-e 1 / 2 +Q(e)) fc _ cj-^fce 1 / 2 +Q(A; 2 e) qj2fc 2 e+0(fce)+0(A; 3 e 3 / 2 ) \ 

l-o;(l-2e 1 /2 + o(£))'= ~~ l-w+w2fce 1 /2 y L + l-w+u^ifce 1 / 2 J 

(w-tjfce 1 /2 + o(fc 2 e))(l-tj+w2fce 1 /2 +CJ 2A: 2 e+0(fce)+0(A; 3 e 3 / 2 )) 

Likewise, the second term from equation (|59|) can be similarly estimated and shown to be 



-he 1 / 2 _ (w-ujke 1 / 2 +O(k 2 e))(l-L0+u J 2ke 1 / 2 +uj2k 2 t+O(k 3 e 3 / 2 )) 
l_ we -2fe e 1 /2 — (l-aj+w2fce 1 /2)2 • 



Taking the difference of these two terms, and noting the cancellation of a number of the terms in 
the numerator, gives (j60j) . 



To see that the error in (I60p is bounded after the summation over k in the range {— e 1//2 , . . . , e 1 ^ 2 }, 
note that this gives 

1/2 ^e 1 / 2 (2feV 2 ) 2 (2fc. 1 / 2 ) 3 _ rl (2t) 2 (2t) 3 , 

Z^fc=-e- 1 /2 l_ a;+tJ (2fc e l/2) (l-w+a;(2fce 1 /2))2 J-l l_ u + w 2i (l-w+w2t) 2 Ut - 

The Riemann sums and integrals are easily shown to be convergent for our uj which lies on f2, which 
is roughly the unit circle, and avoids the point 1 by distance e 1 / 2 . 

Having completed this first step, we now must show that the Riemann sum for the integral in 
equation (j58j) converges to the integral. This uses the following estimate, 

Ee " 1/2 .«.-*« 

fc=-e-l/2 £ max (fc-l/2)£ 1 /2< i <(fc + i/2)e 1 /2 



-2fee 1 / 2 1— 0)e 



< C (61) 



1— we 

To show this, observe that for t G e 1/2 [fc - 1/2, + 1/2] we can expand the second fraction as 

CJe -^ 1/2 (l + Q( £ l/2)) 
l-u,e~ 2kcl/2 (l-2le 1 / 2 +0{e)) 

where I G [—1/2, 1/2]. Factoring the denominator as 



(62) 



(i-^- 2 ^ 2 )( i + -- 2fc ; 1/2 ^t Q( - )) ) (63) 



we can use 1/(1 + z) = 1 — z + 0(z 2 ) (valid since |1 — we 2fcel/2 | > e 1 / 2 and \l\ < 1) to rewrite (|62p 

as 

ce-- 1/2 (l + (e V 2)) (l - ^^gg^ ) 
1 — we ZK<E 

Canceling terms in this expression with the terms in the first part of equation (|6ip we find that we 
are left with terms bounded by 



0{e x / 2 ) Oie 1 / 2 ) 
l_ UJe -2fc e 1 /2 (1 _ we -2fc e l/2 ) 



These must be summed over k and multiplied by the prefactor e 1 / 2 . Summing over k we find that 
these are approximated by the integrals 

£l/2 f-l l-J+u2t dt > f-1 (l-uj+uj2t) 2 dt 

where |1 — u\ > e 1 / 2 . The first integral has a logarithmic singularity at t = which gives | log(l— u)\ 
which is clearly bounded by a constant time | log e 1 / 2 ) for wel!, When multiplied by e 1 / 2 this term 
is clearly bounded in e. Likewise, the second integral diverges like |1/(1 — u))\ which is bounded 
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by e 1 I 2 and again multiplying by the e 1//2 factor in front shows that this term is bounded. This 
proves the Riemann sum approximation. 

This estimate completes the proof of the desired bound when z = 1 + e 1 / 2 . The general case of 
\z\ = 1 + e 1 / 2 is proved along a similar line by letting z = 1 + pe 1 / 2 for p on a suitably defined 
contour such that z lies on the circle of radius 1 + e 1 / 2 . The prefactor is no longer e 1 / 2 but rather 
now pe 1 / 2 and all estimates must take into account p. However, going through this carefully one 
finds that the same sort of estimates as above hold and hence the theorem is proved in general. □ 

This lemma completes the proof of Proposition [27] □ 

Proof of Proposition \28[ We will focus on the growth of the absolute value of the determinant. 
Recall that if K is trace class then | det(/ + K)\ < exp Furthermore, if K can be factored 

into the product K = AB where A and B are Hilbert-Schmidt, then ||-fC||i < ||A||2||-B||2- We will 
demonstrate such a factorization and follow this approach to control the size of the determinant. 

Define A : L 2 (f c ) -> L 2 {f v ) and B : L 2 (T V ) -> L 2 (T C ) via the kernels A((,fj) = e "^°' and 

B(fj ?) _ c IIm(ni c -f(C a -fj3)+ai 2 l/ 3 7 r(-A) z ' V^^^XT^) 
D \'h^) — e e z sin(vr5) r(2 1 / 3 ^'-2 1 /3xT- 1 /3)' 

where z = 2 1/ ' 3 ((' — fj). Notice that we have put the factor e - ' 1 ™'-''^ into the A kernel and removed 
it from the B contour. The point of this is to help control the A kernel, without significantly 
impacting the norm of the B kernel. 



Consider first ||^4||2 which is given by 



A\\ 2 = I I d(dfj-~ 



c 



-2|Im(C)| 



\C-fj\ 2 

The integral in fj converges and is independent of ( (recall that |£ — fj\ is bounded away from zero) 
while the remaining integral in £ is clearly convergent (it is exponentially small as C goes away from 
zero along T^. Thus ||A||2 < c with no dependence on p at all. 

We now turn to computing | \B\ I2. First consider the cubic term £ 3 . The contour is parametrized 
by — + C3«r for r S (—00,00), that is, a straight up and down line just to the left of the y axis. 
By plugging this parametrization in and cubing it, we see that, Re(£ 3 ) behaves like |Im(£)| 2 . This 
is crucial; even though our contours are parallel and only differ horizontally by a small distance, 
their relative locations lead to very different behavior for the real part of their cube. For fj on the 
right of the y axis, the real part still grows quadratically, however with a negative sign. This is 
important because this implies that \e~~s ^ ~ f> )| behaves like the exponential of the real part of the 
argument, which is to say, like 

e -f(|Im(C)| 2 +|Im©| 2 )_ 

Turning to the p term, observe that 

= e Re[(l°g|Al+«arg(-/i))(-Re(5)-iIm(S))] — g - log |/i|Re(5)+arg(-/i)Im(5) _ 
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The esc term behaves, for large Im(z), like e -7r ' Im ( 2 -" . Finally, we must show that the Gamma 
functions have sub-quadratic growth on vertical lines. This follows from Colollary 1.4.4 of [3] which 
states that for x = a + ib and a\ < a < 02, as |6| — > 00 

\T(a + ib)\ = V2^\b\ a - 1/2 e- wW/2 (l + 0(1/|6|)). (64) 

Putting all these estimates together gives that for £ and fj far from the origin on their respective 
contours, \B(f]X)\ behaves like the following product of exponentials: 

e |Im(Ol e -!(|Im(0| 2 + |Im(^)| 2 ) e -^ 

Now observe that, due to the location of the contours, — Re(z) is constant and less than one (in fact 
equal to 1/2 by our choice of contours). Therefore we may factor out the term e~ log ^^ e ^ = \{l\ a 
for a = 1/2 < 1. 

The Hilbert-Schmidt norm of what remains is clearly finite and independent of p,. This is just due 
to the strong exponential decay from the quadratic terms — Im(£) 2 and — Im(r/) 2 in the exponential. 
Therefore we find that ||-B||2 < c l/^| a f° r some constant c. 

This shows that ||-K"a SC,r ||i behaves like \ jl\ a for a < 1. Using the bound that | det(I + Ka SC ' r )\ < 
e \\ K a II we fj n( j th a t 1 det(I + Ka SC,r )\ < e'^'". Comparing this to we have our desired result. 
Note that the proof also shows that K^ C ' T is trace class. □ 

3.4.3. Proof from Section\EM 

Proof of Lemma\3SX We may expand ^(rj) into powers of e with the expression for r/ in terms of 
k{6) from (1511) with a = —1/2 (similarly 1/2 for the £ expansion). Doing this we find 

Re(* ( v )) = e- 1 ' 2 (-le-^Ta cot 2 (%) + ±T [a + k(9)] 2 cot 2 (%)) +0(1). (65) 

We must show that everything in the parenthesis above is bounded below by a positive constant 
times 1 77 — £| for all r\ which start at roughly angle le 1 ^ 2 . Equivalently we can show that the terms 
in the parenthesis behave bounded below by a positive constant times \tt — 0\, where 9 is the polar 
angle of rj. 

The second part of this expression is clearly positive regardless of the value of a. What this suggests 
is that we must show (in order to also be able to deal with a = 1/2 corresponding to the £ estimate) 
that for 77 starting at angle le 1 / 2 and going to zero, the first term dominates (if I is large enough). 

To see this we first note that since a = —1/2, the first term is clearly positive and dominates for 6 
bounded away from ir. This proves the inequality for any range of r/ with 6 bounded from ir. Now 
note that for 6 near n, 

cot 2 (|)«i(vr-0) 2 , tan 2 (f)«4(7r-fl)- 2 , log 2 ( <T= ^ m ) » ^(vr - 6)\ 

We may expand the square in the second term in (|65p and use the above expressions to find that 
for some suitable constant C > (which depends on X and T only), we have 

Re(tf fa)) = e- 1 ' 2 (-i e -V2 Ta (7r - Of + C(ir - 6) 2 ) + 0(1). 
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Since, for some constant c, c — rj\ < (tt — 9) < c|£ — rj\, the second part of the Lemma follows 
from the above equation. Now use the fact that tt — 9 > Ze 1 / 2 to give 

Re(*(r/)) = e- 1 ' 2 {-±lTa(n - 0) + £(tt - 9) 2 ) + O(l). (66) 

Since tt — 9 is bounded by tt, we see that taking I large enough, the first term always dominates 
for the entire range of 9 G [0, tt — Ze 1 / 2 ]. Therefore since a = —1/2, we find that we have have the 
desired lower bound in e" 1 / 2 and \tt — 9\ for the first part of the Lemma. 

Turn now to the bound for Re( x I / (C))- In the case of the rj contour we took a = —1/2, however 
since we now are dealing with the £ contour we must take a = 1/2. This change in the sign of 
a and the argument above shows that equation ([66]) implies the desired bound for Re( x I / (C)), for Z 
large enough. □ 



Before proving Lemma [33] we record the following key lemma on the meromorphic extension of 
z). Recall that z) has poles at \i = t j , j G Z. 

Lemma 45. For fi ^ r 3 , j G Z, nf(fi,z) is analytic in z for 1 < \z\ < t _1 and extends analytically 
to all z ^ or r k for k G Z. This extension is given by first writing z) = g+(z) +g~ (z) where 

f— ' 1 — r K a f— ' 1 — T ix 

k=0 r k=l r 

and where g+ is now defined for \z\ < r _1 and g_ is defined for \z\ > 1. These functions satisfy 
the following two functional equations which imply the analytic continuation: 

9+{z) = — - — + fj.g + (rz), g-(z) = — h -g-(z/r). 

1 — TZ 1 — Z fl 

By repeating this functional equation we find that 

N k N-l _ k 

= Etzt^ + ^ N 9 + (r N zi g _(z) = Y^r z +m~V(^). 

k=l k=0 



Proof. We prove the g+ functional equation, since the g- one follows similarly. Observe that 

oo _. 

g + {z) = ^( T z) fc (l + - 



LIT" 

k=0 r 

OO o t. 



+ ^ i * ( rz ) fc = T - ^ — + M+(tz), 

z — ' 1 — ilT" 1 — TZ 



1 — TZ ^ 1 — LIT 

k=0 ^ 

which is the desired relation. □ 



Proof of Lemma EH Recall that ft lies on a compact subset of C and hence that 1 1 — fiT k | stays 
bounded from below as k varies. Also observe that due to our choices of contours for rj and £, 
\C/v'\ stays bounded in (l,r _1 ). Write z = C/v'- Split jj,f((l,z) as g+(z) + g~(z) (see Lemma Hoi 
above), we see that g+ (z) is bounded by a constant time 1/(1 — tz) and likewise g~(z) is bounded 
by a constant time 1/(1 — z). Writing this in terms of £ and 7/ again we have our desired upper 
bound. □ 
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Proof of Lemma\34\ Observe that uq ~ ^e _1 / 2 loge and hence due to the choice of the contours 
for rj and £, Re(log(£/?7)) = 0(e 1 ^ 2 ). This implies that this first term diverges in e like loge which 
is clearly beaten by e" 1 / 2 (which is a clear lower bound for the right hand side for the choice of 
contours). 

Now consider the infinite sum of logarithms. The closer that l+r n ( gets to zero, the worse the sum, 
so let us assume that the denominator is smaller than the numerator. Then due to the restriction 
on where rj and £ lie on their respective contours we are assured of having |(1 + r n ?])/(l + r n £)| 
bounded above by a constant times e~ 1 / 2 r n |C — tj\. This constant goes to zero as lo increases. 
Summing in n then gives an upper bound of c\( — r/|e -1 and the triangle inequality completes the 
proof. □ 



Proof of Lemma 37, The discussion from equation (143D to (I44j) is rigorous, and, for rj and ( in the 
fixed compact regions of their contours, we have uniform control over the error in e and uniform 
convergence of the q-Gamma functions to standard Gamma functions. This suffices to prove the 
Lemma. □ 



Proof of Lemma 



Recall that we have defined 

- 1 '\-a + 



m 



The argument now amounts to a Taylor series expansion with control over the remainder term. Let 
us start by recording the first four derivatives of A(£): 



HO 
A'(C) 

A"(C) 
A"'(C) 
A""(C) 



X 

i-C 



xlog(l-C) + +mlogC 
+ 



(i-C) 



t 1 in 



C 



X 


+ 


2t 




(i-C) 3 


2x 


+ 


m 


~ (i-C) 3 


(i-C) 4 


Gx 


+ 


24t 


~ (i-C) 4 


(i-C) 5 



+ 



m 

2m 

6m 



We Taylor expand ^(C) = A(£) — A(£) around £ and then expand in e as e goes to zero and find 
that 

| e -l/2 + (l) 

f £ - 3/2 + o(r') 

0(e" 3 / 2 ). 



A"(£) 
A"'(£) 
A""(£) 



A Taylor series remainder estimate shows then that 

l*(c) - [a'(o(c - + aa w (o(c - e) 2 + aa w (o(c - 6 3 ] I < sup i|A""(t)nc - £i 4 , 

t6B(e,ic-€i) 

where -B(£, |C — £|) denotes the ball around £ of radius |£ — £|. Now considering the scaling we have 
that £ — £ = c^" 1 e 1 / 2 C so that when we plug this in along with the estimates on derivatives of A at 
£, we find that the equation above becomes 



*(C) 



,1/3 



Ole 1 / 2 ). 
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The above estimate therefore proves the desired first claimed result. 

The second result follows readily from \e z — e w \ < \z — w\ max{|e 2 |, \ e w \} and the first result, as well 
as the boundedness of the limiting integrand. □ 



Proof of Lemma \39[ Expanding in e we have that 

. = ^^ = l-^» + 0(.) 

where the error is uniform for our range of fj' and ( and where z = C3 1 (( — fj'). We now appeal to 
the functional equation for /, explained in Lemma [4"5l We wish to study e 1 > 2 g+(z) and e 1 ^ 2 g-(z) 
as e \ and show that they converge uniformly to suitable integrals. First consider the g+ case. 
Let us, for the moment, assume that < 1. We know that \rz\ < 1, thus for any N > 0, we have 

,l/2„ ( \ _ 1/2 V^AT fi k 1/2 r.N „ (~N ~\ 

Since, by assumption, < 1, the first sum is the partial sum of a convergent series. Each term 
may be expanded in e. Noting that 

1 _ T k z = 1 - (1 - 2e 1/2 + 0(e))(l - e 1/2 z + O(e)) = (2k + z)e 1/2 + fcO(e), 

we find that 

6 1- T k z — 2k+z ^ ftlV V e )• 

The last part of the expression for g + is bounded in e, thus we end up with the following asymptotics 

e 1/2 9 + (z) = E£Li A + ^O^ 2 ) + AMI)- 

It is possible to choose iV(e) which goes to infinity, such that A^ 2 0(e 1 / 2 ) = o(l). Then for any 
fixed compact set contained in C \ {—2, —4, —6, . . .} we have uniform convergence of this sequence 
of analytic functions to some function, which is necessarily analytic and equals 

Eoo fi k 
fc=l 2k+z- 

This expansion is valid for < 1 and for all z G C \ {—2, —4, —6, . . .}. 

Likewise for e 1 / 2 g-(z), for \p,\ > 1 and for z £ C \ {—2, —4, —6, . . .}, we have uniform convergence 
to the analytic function 

-,k 



E 



fi K 



-oo 2k+z 



We now introduce the Hurwitz Lerch transcendental function and relate some basic properties of 
it which can be found in |37j . 

for w > real and either \a\ < 1 and s £ C or \a\ = 1 and Re(s) > 1. For Re(s) > it is possible 
to analytically extend this function using the integral formula 

where additionally a £ C \ [l,oo) and Re(w) > 0. 



44 



I. CORWIN AND J. QUASTEL 



We can express our series in terms of this function as 

£?=i 2& = 1, 1 + 5/2), EL-oo #1 = -^(A- 1 , 1, -z/2). 

These two functions can be analytically continued using the integral formula onto the same region 
where Re(l + z/2) > and Re(— 5/2) > - i.e. where Re(z/2) G (-1,0). Additionally the analytic 
continuation is valid for all jl not along M + . 

We wish now to use Vitali's convergence theorem to conclude that JXf(JX,z) converges uniformly 
for general jl to the sum of these two analytic continuations. In order to do that we need apriori 
boundedness of e 1//2 <7+ and e 1//2 g_ for compact regions of jl away from M + . This can be shown 
directly as follows. By assumption on jl we have that |1 — r k jl\ > c _1 for some positive constant c. 
Consider e 1 / 2 ^ first. 

We know that |rz| is bounded to order e 1 / 2 away from 1 and therefore this shows that |e 1 / 2 5+(z)| 
has an upper bound uniform in jl. We can do a similar computation for e 1 / 2 g-(z) and find the 
same result, this time using that \z\ is bounded to order e 1//2 away from 1. 

As a result of this apriori boundedness uniform in jl, we have that for compact sets of jl away from 
M + , uniformly in e, e 1//2 g+ and e l l 2 g- are uniformly bounded as e goes to zero. Therefore Vitali's 
convergence theorem implies that they converge uniformly to their analytic continuation. 



Now observe that 



and 



^,l,l + z/2) = Uo° ^dt 



-i$(/r\ 1, -2/2) = -U °° "-^Wdt = I 
Therefore we can combine these as a single integral 

1 roo £e"~ a/2 A4- — 1 r°° fis'*/ 2 ds 

2 J — oo e'— /t 2 JO s— /i s ' 

The first of the above equations proves the lemma, and for an alternative expression we use the 
second of the integrals (which followed from the change of variables & = s) and thus, on the region 
Re(z/2) 6 (—1,0) this integral converges and equals 

\ir(-jl)~ z csc(7rz/2) 

which is analytic for jl E C \ [0, oo) and for all z E C \ 2Z. Therefore it is the analytic continuation 
of our asymptotic series. □ 



4. KPZ EQUATION LIMIT OF WASEP 



Proof of Theorem^ First let us describe in simple terms the dynamics in T of Z e (T, X) defined in 
(|13p . To ease the presentation, we will now drop all subscripts e. There is a deterministic part and 
random jumps. The jumps are at rates 

r_(x) = e~ 2 g(l — rj(x))r)(x + 1) = -e~ 2 q(l — fj(x))(l + rj(x + 1)) 
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to e~ 2X Z and 



r + (x) = e 2 pr/(x)(l — rj(x + 1)) = -e 2 p(l + r)(x))(l — f/(x + 1)) 



to e 2X Z, independently at each site X G eZ. We write it as follows, 

dZ = QZdT + (e~ 2A - \)ZdM- + (e 2A - l)ZdM+ 

where 

Q = e~ 2 v + (e" 2A - l)r_ + (e 2A - l)r+ 

and dM±{T,X) = dP±{T,X) - r±{X)dT where P_(T, X), P + (T, X), I e fZ are independent 
Poisson processes running at rates r_(T, X),r+(T, X). Let 

D = 2^=l-^e + 0(e 2 ) 

and A = A £ be the eZ Laplacian, A/(X) = e- 2 (f(X + e) - 2f(X) + f(X - e)). We also have 

§DAZ = ±e- 2 D(e- Xf i( x+ V _ 2 + e A ^)Z 

The parameters have been carefully chosen so that 

n = ie- 2 D{e- Xflix+1) - 2 + e x ^ x) ) 

We can do this because the four cases, corresponding to the four possibilities for fj(x),fj(x + 1): 11, 
(— 1)(— 1), 1(— 1), (—1)1, give four equations in three unknowns, 

11 ±e~ 2 D(e~ A -2 + e A ) = e~ 2 u 

(-1)(-1) le- 2 D(e X -2 + e- x ) = e- 2 v 

1(-1) \(- 2 D{e x - 2 + e A ) = e" 2 ^ + (e 2A - l)e~ 2 p 

(-1)1 \e- 2 D{e- x - 2 + e~ A ) = e~ 2 ^ + (e~ 2A - l)e- 2 q 

But luckily, the first two equations are the same, so it is actually three equations in three unknowns. 
This is the Gartner transformation [20]. The solution is A = \ log(q/p), D = 2^/pq, v = p+q—2^/pq. 



Hence [2D!,P2], 

dZ e = ±D e A e Z e dT + Z e dM e (67) 

where 

dM e (X) = (e~ 2Ae - l)dM-(X) + (e 2Ae - l)dM+(X) 
are martingales in T with 

d{M e (X),M e (Y)) = e~ l l{X = Y)b e (T_ [t -i x]V )dT 
where T x r](y) = r](y — x) and 

b e { V ) = l-f)(l)fj(0) + b £ (r l ) 

where 

Uv) = e- 1 {[p((e- 2AE -l) 2 -4e) + (7 ((e 2Ae -l) 2 -4e)] 
+ [q{e- 2X ^ - I) 2 -p{e 2X < - 1) 2 ](^(1) - r?(0)) 
-[q(e- 2X ^ - l) 2 +p(e 2X * - l) 2 - e]^(l)^(0)}. 
Clearly b t > 0. It is easy to check that there is a C < oo such that 

< L < Ce 1/2 



16 



I. CORWIN AND J. QUASTEL 



and, for sufficiently small e > 0, 

< b e < 3. 

We have the following bound on the initial data. For each p = 1,2,... there exists C p < oo such 
that for all X E R, 

E[Z*(0,X)] < e CpX . (69) 

For any 5 > let denote the distribution of Z e (T, X), T G [5, oo), as measure on D[[5, oo), C(R)] 
where D means the right continuous paths with left limits, with the topology of uniform convergence 
on compact sets. In [12], Section 4 it is shown that if (|69l) holds, then, for any 5 > 0, e > 
are tight. The limit points are consistent as 5 \ 0, and from the integral version of (167|) . 



Z t (T,X) = e Pe(T,X-Y)Z e (0,Y) 

+ f e Pe(T - S,X -Y)Z e (S,Y)dM e (S,Y) 



Y£e2 

where p e (T,X) are the transition probabilities for the continuous time random walk with generator 
±DA t , normalized so that p e (T,X) — > p(T,X) = e /^f ' we ^ ave 



lim lim E 



(z e (T,X)-eJ2 Pe(T,X -Y)Z t (0,Y)\ 

V Y£tZ J 



so that the initial data (1151) hold under the limit . Finally, we need to identify the limit of the 
martingale terms. Recall the key estimate in [12] which, translated to our context, says that for 
any < 5 < Tq < oo and p > 0, there are C\, C*2 > such that for all 5 < S < T < Tq and e > 0, 

E[\E[(Z e (T,X + e) -Z e (T,X))(Z e (T,X) - Z e (T, X - e))\^{S)}\] < e 1/2 " p |r - 5|- 1/2 e a|x| . (70) 



Again, this is only proved using (|69|) . and extends without change to the present context. Let us 
briefly recall why such a thing is true. It is well-known in the theory of stochastic partial differential 
equations that the solutions of a stochastic heat equation will be Holder 1/2 — p in space, for any 
p > 0. This is proved using only the integral version of the equation and L p bounds on the initial 
data. Hence it extends in a standard way to a discretization such as (|67p of such an equation, as 
long as we have (|69p . with constants independent of e. 

Let <*p be a smooth test function on R. We hope to show that under 

N T ((f):= [ ip(X)Z(T,X)dX - | f [ (p"(X)Z(S,X)dXdS 
Jr Jo Jr 

and 

A T (tp) := N T {^f - \ [ l\ 2 {X)Z 2 {S,X)dXdS (71) 
Jo Jr 

are local martingales. Because we also have E[Z 2 (T, X)] < e c \ x \ for all T > 0, we have uniqueness 
of the corresponding martingale problem, following Section 5 of [T2] and Theorem 2.2 of [TU]. That 
Nx is a local martingale under & follows from (|67p which says that microscopically, 

N T , e (<p) ■= [ y(X)Z e {T,X)dX-\ [ [ D t A t (p(X)Ze(S,X)dXdS 
Jr Jo Jr 
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is a martingale under The key point is to identify the quadratic variation, ie. the martingale 
Ay ((/?). Microscopically we have that 



in probability. Together with (j68|) this shows that At (9?) defined in ([7T]) is a martingale under 



5. Manipulations and asymptotics of the edge crossover distributions 
We now consider various asymptotics and properties of the edge crossover distributions. 



5.1. Large T asymptotics of the edge crossover distributions (Proof of Corollary [5]). 

The proof of Corollary [5] proceeds similarly to that of Corollary 3 of [2 J . The first step is to cut 
the /2 contour off outside of a compact region around the origin. Proposition 18 of [2] (with the 
modifications explained in Section I3.2f) shows that for a fixed T, the tail of the \x integrand is 
exponentially decaying in p,, and a quick inspection of the proof shows that increasing T just serves 
to speed up this decay. This shows that we can cut off the infinite tails of the C contour at cost 
which goes to zero as the cut occurs further out. 

From now on we may assume that /2 lies on a compact region along C. If we plug in the scalings 
for space as 2 1//3 T 2//3 X and fluctuations as 2~ 1//3 T 1 / 3 s and make the change of variables that 
£ = T~ 1//3 £,77 = T^ 1 / 3 ^, and fj = T" 1 / 3 ^' then we find that the integrand in the kernel for K^ ge 
can be written as 



can be deformed back to their original form by a combination of Cauchy's theorem and Proposition 
1 of |H], which says that as long as we don't pass any poles in the kernel, we can deform the contours 
on which an operator acts without changing the value of the Fredholm determinant. 

Let and T v be these rescaled contours. The only requirement on these contours is that they 
look like those given in Figure [2] and that they both go to the right of the right most pole of the 
Gamma functions which occurs at X. 




Following the argument in Section 4 of |12j . shows that (|70p suffices to prove that 




which completes the proof. 



□ 
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We now claim that with the scalings described above det(I — Ka ds °) L 2^ \ converges, uniformly in 
fb, to det(I - K a ) L 2( Fv ) where K a has kernel K a (r],r]') given by 

jy~ / / \ f exp{-|(C 3 - rf 3 ) + s(( - rf) + X 2 (( - rf)} rf-X 

K ' M = k (C-MC-rf) C^" C (?3) 

The claim follows exactly as in the proof of Corollary 3 of [2] and relies on the fact that in the 
scalings present in (|72p the first fraction converges for compact sets of £ and rf to 1/(C — v')? while 
the second fraction converges to (7/ — X)/(<^ — X). Convergence for compact sets of Q and rf is 
enough since the exponentials provide sufficient decay for the necessary trace class convergence of 
operators. 

Now we use the method of [H] to factor (|73p into a product of three operators ABC and then 
reorder as BCA without changing the value of the determinant. Observe that given our choice of 
contours, Re(£ — i]') < and hence 

exp{s(C - r/)} 



/oo 
exp{x(( — r}')}dx. 



C-77' 

Inserting this into fl73]) we find that K a = ABC where A : L 2 (s,oo) -> L 2 (T^), B : L 2 (T C ) -> 
L 2 (s,oo) and C : L 2 (F V ) — > L 2 (F^) and are given by their kernels 

A( V ,x) = exp{±r/ 3 -(x + X 2 )r/}(? ? -X) (74) 

B(x,Q = e W {lC 3 -(x + X 2 )C} (75) 

CM = (C-vK-xy (76) 

Reordering does not change the value of the determinant, and we are left with an operator 
BCA(x,y) acting on L 2 (s, 00) with kernel 

BCA*,y)—f f "P<-C 3 /3 + gV3 + ye - g ± XH< - „)} n _ - X 
Jt c Jt v C-v C-x 

Shifting the x and y contours by X 2 and using = 1 + f^Y- we nave an operator BCA acting 
on L2(s — X 2 , 00) with kernel 

una \ I I ™P{-C 3 /3 + .) 3 /3 + K-"rf (,. 2-C \ AA . 

bcm,.v) — J r( j Ti — v (i + r ^)m. 

Expanding the multiplication, the first term corresponds to the Airy2 kernel, 

/•oo 

K A2 (x,y)= / Ai(t + x)Ai(t + y)dt. 
Jo 

The second term is 

exp{-C 3 /3 + 77 3 /3 + yC-a^} , 
^ 

We can factor this into the rj and £ integrals separately. The r\ integral gives Ai(x). The £ integral 
can be evaluated as follows. First recall that due to the dimple in the definition of IV, £ — X is on 
a contour which lies to the right of the origin. Make the change variables to let Z = £ — X, which 
gives for the £ integral, 

X 3/3+Xy 



J d Ze- z3 / 3 - bz2+cZ ^, 
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where b = X and c = —X 2 + y. We now appeal to Lemma 31 (B) of Appendix (A) of [6] which 
states that (recall we have absorbed factors of (2iri)~ 1 into our dZ) 

J dZe- z3 /*- bz2+cZ ^ = -e~l bA - bc J™ dtAi(b 2 + c + t)e~ bt . (77) 

The integral in the above equation, however, is over a contour which is to the left of the origin, 
where as our integral is to the right. This is easily fixed by deforming through the pole at Z = 
which gives 

-X 3 /3+Xy fi,f jy^-Z^/Z-bZl+cZ 1 



e \ 1 + J dZe ~ Z ) ' 

where the Z integral is now to the left of the origin. Applying (|77p we are left with 

(roc \ /-oo 

1 - e * 3 /3-A-, J dtM ( y + t)e ~Xt\ = e -XV3 + Xy _ j dtM{y + t)e -Xt 

Thus the final kernel is 

(/•OO 
e -X*/3+Xy _ J dtM{y + t)e -Xt 

which one can recognize from Definition [20] of Section 12.21 as 

5.2. Alternative forms for the edge crossover distributions. In this section we develop an 
alternative formula for the edge crossover distribution. Our starting point is equation (|48p for 
A"f c ' r . This can be transformed into Kf gc by taking a = 2- 1 / 3 T l / 3 s and X = 2 1 / 3 T 2 / 3 X'. We 
will stick to the original form however. Recalling the equation, we have a kernel 



1 /v3 ~/3^ , nl/3-.r? ~l\ I ol/3 / / ^ e jj. \ 1 \ Z _± T> d( = 

r(2Wfr-£)c-ij 



exp |-|(C 3 - rf ) + 2V3 o(C _ jfl j 2 y^J ^ ^ - _ - dtj 
For Re(z) < we have: 



oo 

xz 



e xz dx = 

z 



which, noting that Re(£ — fj) < 0, we may apply to the above kernel to get 

-2 2 / 3 / r fexp {-^ ( C 3 - t) ~ 2^aA M^^l ^ ~ f ) 

This kernel can be factored as a product ABC where 

A : L 2 (a, oo) -> L 2 (f^), 5 : L 2 (f ? ) -> L 2 (a, oo), C : L 2 (f „) L 2 (f c ), (78) 
and the operators are given by their kernels 

A(fj, x) = e 2 1/3 (-^, B (x, C) = e 2l/3 ^, 

x 9rt /"°° f ^ 1 /ie- 2173 ^') r(2V3£ - &) 



T - 

2/ 



Since det(/ — ABC) = det(I — BCA) we consider BCA acting on L (a,oo) with kernel 

f 00 Jn2/3 /" /" „„„ f T /a3 ~3\ , ol/3/„ , ol/3/„, , r (2 1 / 3 C - ^) 



e 1 — jl 



22/3 /r- /r- 6XP " ^ + 2l/3(X + * )C ~ " 2V3(y + m ] r(2i4-|) d ^ C " 
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Recall the two integral formulas: 

poo 

T(z) = / s{- l e- sl d Sl (79) 
Jo 

where C is counterclockwise from oo to oo going around R + (the branch of the logarithm is cut along 
K+). Both equations holds for Re(z) > and can be analytically extended using the functional 
equation for the Gamma function. We can rewrite 



L L 6XP {"l (C ~ 3 " f) + 2V3(X + t)C " " 2V3(2/ + m } r^4_|) ^ 
itiJcT i J r ex p{-|(C 3 -0 3 ) + 2 1/3 (x + t + log Sl )C-2 1 /3( 2/ + t + i O g(_ S2 ))^| 



— 2£. x 

xdfjdQs 1 T {-s 2 )~t e~ Sl ~ S2 dsid-S2 



Using the formula for the Airy function given by 

Ai(r) = I exp{— -z 3 + rz}dz 

we find that our kernel equals 

_ l 2 2/3 T -2/3 f r r Ai(r- 1 /3 2 V3 (a . + i + logSl) ) Ai ( r -l/3 2 l/3 ( +t + 1 ( _ S2)))rfi 



^ T" 1 



xs x T (-s 2 ) T e Sl S2 dsids 2 . 

Note that the formula is only making sense as an integral for X > and has to be extended by 
analytic continuation to other values. 

We can also write it in compact form as follows. Define the Gamma transformed Airy function 

Ai r (a, b,c)= [ exp{--z 3 + az}T(bz + c)dz (81) 



- c 

and the inverse Gamma transformed Airy function 



f 1 1 

Ai r (a,6,c)=/ exp{-z 3 -az} — — — dz 
Jt v 3 T(bz + c) 



(82) 



These are automatically well defined, as long the contours are such that for a, b, c they avoid the 
poles of the Gamma function. Plugging these in we get 

jf jf exp |-|(C 3 - rf ) + 2 1 /3( X + t)c - 2 V3 (y + E|^|z|| d ^ ( 83 ) 

= T^Ai^/e^s + t), k^ 1 , --)Ai r (^ 1 (2/ + >^ ~f ) (8 4 ) 
where k t = 2~ 1 / 3 T 1 / 3 . Thus we find 

/OO ~ j£ 

_ t _ - A.i r (K^{x + t), k^ 1 , -— )Ai r (^ 1 (y + t), k^ 1 ,-—). (85) 
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5.3. Airy Gamma asymptotics. In this section we obtain some asymptotics of the Gamma 
transformed Airy functions (|8ip and (|82p which will prove useful in Section 15.41 We start by 
observing some bounds on Gamma functions. Recall the functional equation 

T(z + 1) = zT(z) (86) 

and the bound (0], 6.1.26) 

\T(x + iy)\ < \T(x)\. (87) 

The following two asymptotic bounds are standard (see [3] 1.4.3 and 1.4.4 respectively). For 5 > 
and | arg z\ < ir — 6, 

T(z) = V2^z z - 1/2 e~ z (l + o(l)) as|z|^oo. 



When z = x + iy and x\ < x < xi and \y\ — > oo, then 

\T(z)\ = V2^\y\ x - 1/2 e~^ (1 + 0(l/|y|)) , (89) 
where the constant implied by 0(l/|y|) depends only on x\ and xi- 

Lemma 46. There exists a constant C > such that for allr > and all z = e ld for G [— 37r/4, 37r/4] 
we have 

T(rz) < CT(r). (90) 

Proof. Consider separately the three cases: (1) when r < 1/2, (2) when r G (1/2,2) and (3) and 
when r > 2. For case (1), when r < 1/2 we may apply the functional equation (|86p once giving 
T(z) = T(z + l)/z. Since \z\ < 1/2, it is clear that \z + 1| > 1/2. Therefore, using (|57|) we find that 

|r(*)| < ^ < g Of) 

where x is real and bounded in [1/2,3/2], and hence the bound T(x) < C for C = T(3/2). Since 
r(r) ~ c/r for r < 1/2 it follows that we can express the bound C'/\z\ in terms of CT(r) for an 
appropriate constant C - as desired. 

Turning to the case (2), when r e (1/2,2) we may apply the functional equation k times so that 
1/2 < Re(z + A;) < r. Here k is at least 1, but at most 3. This shows that 

Since arg(z) G (— 3-7r/4, 37r/4) it follows that \z + k\ > c for a fixed positive constant. Thus, taking 
absolute values in the above equation and using this bound we find 

\T(z)\ < C\T(z + k)\ < CT(r) (93) 

where C is a fixed constant (bounded by 1/c 3 ). 

Turning to case (3), when r > 2 we may apply the functional equation k times so that r — 1 < 
Re(z + k) < r. Here k > 3 and since arg(z) G (— 37r/4, 37r/4) it follows that \z + k\ > 1. Thus 
taking absolute values in (|92p we have 

\T(z)\ < \T(z + k)\ <T(r), (94) 

as desired. □ 
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Here is another, rather weak inequality we will use. 

Lemma 47. There exists a constant C > such that for all r > we have 

|r(re 3i " /4 )| < C/r (95) 

Proof. For r < 2 this follows immediately from Lemma [4U1 For r > 2 set z = re 3 " 1 / 4 . We may 
apply the functional equation (|86p k times so that —1/2 < Re(z + k) < 1/2. Here k > 3 and since 
arg(z) = 37r/4 it follows that \z + fe| > 1. Thus taking absolute values in ([92]) we have 

|r(js)| < |r(* + A;)|. (96) 

We may now apply Lemma [89] with x\ = —1/2 and %2 = 1/2. This clearly implies the desired 
decay (actually much stronger decay than needed). □ 

Lemma 48. For all constants c > it/2, there exist C > such that for all z with Re(z) > 

\l/T{z)\<Ce c W. (97) 

Proof. This is established in two parts. For z such that |Im(z)| < 1 this follows from functional 
equation ([So]) and the boundedness of 1/T(z) for < Re(z) < 1 and |Im(^)| < 1. Similarly, for z 
such that |Im(z)| > 1 we may first consider such z which also satisfy < Re(z) < 1. By (]89p . these 
z satisfy (i97|) . This bound can be then be extended to all z with |Im(z)| > 1 and Re(z) > by the 
functional equation. □ 



Recall the definitions (|8ip and (|82p of the Gamma transformed and inverse Gamma transformed 
Airy functions. In the following lemma constants may change line to line. 

Lemma 49. Fix b > 0, then 

(1) There exists a constant C such that for a > 0, 

|Ai r (a,6,0)| < C((l + |a|)- 1 r(6(l + |a|)- 1 ) + 6- 1 ) . (98) 

(2) For any e > there exists a constant C such that for all a < 0, 

|Ai r (a,6,0)| < C ((1 + |a|)- 1 r(6(l + la))" 1 ) + — + j + (1 + |a|) e 6^-^ . (99) 

(3) There exists a constant C such that for all a >0, 

|Ai r (a,6,0)| < Ce~l ai '\\ + |a|)" 1/4 . (100) 

(4) There exists a constant C such that for all a < and all c > ir/2, 

|Ai r (o,6,0)| < Ce cb|a|1/2 . (101) 

Proof. The lemma is intended to give decay bounds as \a\ grows. We can split consideration up 
into \a\ > 1 and \a\ < 1. For \a\ < 1 direct inspection of the integrals reveals the above claimed 
bounds. What follows, therefore, deals with the \a\ > 1 bounds. We will establish these bounds in 
terms of their dependence on \a\. However, in order to state the results of the lemma for all a, we 
then make the modification of replacing \a\ by 1 + \a\ which, up to constants, does not affect the 
validity of the bounds. 
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To prove equation (198j) change variables z i— >■ a 1//2 z and the deform the image of the contour f £ 
so that it is given by three pieces: a ray coming from e~ 3i7r / 4 oo to a positively oriented arc of the 
circle centered at the origin of radius a" 3 / 2 (from angle — 3tt/A to 37r/4) and finally a ray going 
towards e 3l7r/,4 oo. The integral now is 

Ai r (a,6,0) =a 1/2 J exp{-a 3/2 (^ 3 + ~z)}T{ba l/2 ~z)dz. 

Let us first consider the integral along the circular arc. Along this arc, Re(— a 3//2 (|z 3 +z)) is of order 
1. For the Gamma function we may apply Lemma 06] to see that |r(6a 1//2 5)| < CT(&a -1 ). Since 
the arc has length of order a -3 / 2 we find that (recalling the a 1//2 prefactor above) the contribution 
of the circular arc is like Ca~ 1 T(ba~ 1 ). 

We must consider the two rays, though by symmetry it suffices to consider just one. The argument 
of the exponential can be bounded in real part by -a 3 / 2 ~z and the Gamma can be bounded (Lemma 
|47|) by ba i%^\ ■ Thus the integral along the ray may be bounded above by the following real integral: 

3 /2 r Cdr C 



,V2 / e ~ca 
-3/2 



ba l / 2 r b 



Thus the contributions of the integrals along the rays are like C/b, and hence the bound is estab- 
lished. 

To prove equation (|99p replace a by — a and change variables z \— > a}/ 2 z, giving 

Ai r (a,6,0) =a 1/2 J exp{-a 3 / 2 (^ 3 + z)}T{ba}' 2 z)dz. 

For the contour in the resulting integral choose the following: a ray coming from e _3j7r / 4 oo to —i, 
a line segment from —i to — a _3//2 i, a positively oriented arc of the circle centered at the origin of 
radius a -3 / 2 going from — d~ 3 l 2 i to d~ 3 / 2 i, a line segment from d~ 3 l 2 i to i and a ray from i to 
e 3i7r / 4 oo. 

As before, it is easy to show that the integral on the arc is bounded by Cd~ 1 T(bd~ 1 ) for some 
constant C > 0. The integral on the rays can be easily bounded as well. Just as in the proof of 
Lemma H71 we may establish that along the ray 

where the last inequality is from the fact that \z > 1 along the contour. Thus the integral along 
the ray is bounded by 

% [ exp{-d 3 / 2 (-z 3 + z)}dS. (103) 
b J 3 

This, in turn, can be bounded by the real integral 

s~i poo rit 

- / exp{-cd 3 / 2 r}dr < — — (104) 
b J ba 6 / 2 

for some constant C > 0. Thus the integral along the rays are bounded by C/(bd 3 / 2 ). 
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It remains to control the integral along the line segments. By symmetry it suffices to consider just 
the segment from a~ 3 / 2 i to i. By the triangle inequality we can consider the norm of the integrand, 
and along this contour the exponential term is uniformly of norm 1. Thus the integral along the 
line segment is bounded by 

a 1/2 / \T{bd l/2 ~z)\d~z < b' 1 / \T{z)\dz. (105) 

Ja~ 3 / 2 i Jba~ 1 i 

The integral 

root 

\T(z)\dz = 0(r- e ) (106) 



for any e. To check this fact consider the two cases r < 1 and r > 1. If r < 1 then split the integral 
(|106p into two parts - from ri to i and i to ooi. On the first part observer that |T(z)| < C/\z\ (as 
follows from Lemma[l6]and the small r asymptotics of T(r)). This bounding integral with l/|z| can 
be evaluated to equal log(l/r) which is bounded by Cr~ € for any e (here C depends on e though). 
The remaining integral from i to ooi is easily bounded by a constant by using the asymptotics of 
(f89|) . Thus we find the desired bound claimed in ()106j) . When r > 1, the asymptotics of ([89]) easily 
yields the desired bound. 

Using the estimate of (jlQ6f) we may bound (I105D and conclude the desired bound of C|a| e 6 _1_e . 

To prove equation ()100p change variables z t— >■ a 1 / 2 z and choose the integration contour to be 
the following: a ray from e _i7r//3 oo to 1 and then a ray from 1 to e iw ' 3 oo. Along these contours, 
the reciprocal of the Gamma function is bounded by a constant (this can be seen by applying 
the Stirling's formula asymptotics (i88l) ). and hence can be removed, leaving us with the standard 
asymptotic analysis for the Airy function - thus follows the formula. 

To prove equation (jlOip replace a by —a and change variables z i— > d l / 2 z, giving 
Ai r (a, b, 0) = a 1 ' 2 j exp{a 3 / 2 (^z 3 + z)} — zj^ dz. 

We may choose the integration contour to be the following: a ray from e~ t7T ^oo to —i, a line 
segment from —i to i, and a ray from i to e J7r//4 oo. Let us consider the integral along the line 
segment. We may bound the reciprocal of the Gamma function on the imaginary axis by Lemma 
HHl From this bound and the fact that along this line segment, Re(|z 3 + z) = 0, one easily bounds 

the integral by Ce cbal/2 (note that the a 1 / 2 prefactor can be absorbed into the exponential term 
through the constant c). It remains to establish a similar upper bound on the integral over the two 
rays. By symmetry we can consider only the ray going from i to e J7r / 4 . The bound of Lemma H51 is 
valid as long as Re('w) > and thus we may substitute it into the integrand. This yields a bound 
of (we have absorbed the a 1 / 2 as above) 

C J exp{a 3 / 2 (iz 3 + 2)}e cb&1/2 ' z dz, (107) 

with the integral over the ray in question. Due to the rapid decay of Re(|z 3 + z) along the ray, 
it is now easy to show that this is (just as with the integral along the line segment) bounded by 

Ce cbal 2 where still we are assuming c > tt/2. Combining these bounds yields the claimed formula 
(fTOTp . □ 
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5.4. Upper tail of F^ e . Start with the formula 

1 - i^ d ge (s) = - J_e-^[det(I - k T ,a) ~ detl]. (108) 
where the determinants are evaluated on L 2 (s, oo) and 

K T ,n(x,y)= [ _^f_ - Ai r (x + t, K-\0)Ai r {y + t, k~\0) 

J —oo & r 1 

is obtained by rescaling Kx,T,fr is given in (|85[) . First of all note that 

det(/ - k T ,fi) = det(J - A) where A = U^Kt^U. (109) 
We will use this with 

Uf(x) = (x 4 + l)- 1 / 2 f(x). (110) 

We also make use of the fact that if A = A\A 2 for A\ and A 2 Hilbert- Schmidt, then A is trace-class 
with 

|det(/ + A)-det/| < ||A||ie l|A||l+1 < ||Ai|| 2 p 2 || 2 e l|Al||2||A2||2+1 . (Ill) 

We factor A = A 1 A 2 where A 1 : L 2 (R) ->■ L 2 (s,oo), A 2 : L 2 (s,oo) ->■ L 2 (M) are defined by their 
integral kernels as 

Ai(ar,i) = Ai r (i + t^ 1 ,0)(x 4 + l)- 1 / 2 (t 4 + l)- 1 / 2 , (112) 

= e -^f_ + °)^ 4 + 1 ) 1/2 ( t4 + 1 ) 1/2 - ( 113 ) 

We estimate the Hilbert-Schmidt norms with the aid of Lemma [4*91 In what follows, constants are 
denoted by upper and lower case c and can change value from line to line. For A\ we use the bound 

|Ai r (a,6,0)| < C|6|- 3 /V|a| 2 + l- (114) 
which holds for all a, b from Lemma [J9l We thus have 

/OO fOO 
dt dx(\x + t\ 2 + l)(x 4 + l)-\t 4 + l)~ 1 dxdt <CT. (115) 
-oo J —oo 



-oo J —oo 

For Ao we start with the bound 



e -K T t _ 



<C\fl\(e KTt Al), (116) 



as follows from the fact that the contour on which jl lies is bounded from the positive real axis by 
a uniform distance. Then we use (jlOOp and (jlOip to get 

/■oo /'oo 4 

||A 2 ||| < C|/i| 2 / dx dt(e 2KTt Al)e-3 lx+tli/ (x 4 + l)(t 4 + l) 

J S J — X 

/"OO /* — 3^ 

+C dx dt{e 2KTt A lje- 4 ^ 1 !* 4 -*' 172 ^ 4 + l)(i 4 + 1)) := C\fi\ 2 (h + J 2 ). 

7s J —oo 



Look at / 2 first. We are only interested in s > 1, so —x < —s < 0, and therefore e 2KT * < 1 so we 
have 

fOO r—X /-OO i _ 3 

I 2 = (x 4 + l)dx / dte 2KTt ~ 4K T ^=*(t 4 + 1) < / (x 4 + l)dxe- 2KTX ((x 4 + 1) + CV^e^T . 

is J ~oo Js 
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Assuming T > Tq > 0, we have 

h < Cs s e~ 2KTS 

Now for I\, we write it as I\ = I3 + 1 a where I3 denotes the i integration from — x to —x/2 and I4 
denotes the i integration from —x/2 to 00. Now 

roo p — x/2 



/oo p — x/z 

dx dte 2KT \x 4 + \){t 4 + 1) <csV 

J —X 



and 

COO /"OO 



poo poo 4 . 

I 4 < / dx dte~3\ x+t \ 3/2 ( x 4 + l)(t 4 + 1) < cs 8 e 



„8„-cs 3 / 2 

-x/2 

Note that the c in the last exponent can be computed to be about 4-3 -1 -2~ 3 / 2 but it is not optimal 
so we do not pursue it. We obtain 

P2II2 < C\fl\ 2 s 8 (e- 2KTS + e- cs ' i/2 ). (117) 

In summary, combining (|llip . (|115|) . (|117|) we obtain a bound 

I det(/ - K T ,fi) ~ det I\ < C\~^\s 4 T l / 2 {e~ cTl/3s + er cs3/2 ) e I^ (118) 

as long as C s^T 1 / 2 (e- cTl ' 3 s + e" cs3/2 ) < 1/2. Note that c = \ e^+l^l < 10. From (fTUHD . 
we obtain with a new constant C = 10C, 



|1 - F^ ge (s)| < C S 4 T 1 / 2 (e- cTl/3s + e~ cs3/2 ) (119) 

whenever the right hand side is less than or equal to 5. But since the left hand side is less than or 
equal to 1, the inequality holds for any s. 



5.5. Proof of small T edge crossover distribution behavior. 

Proof of Proposition We have 

„ „ n n 

E[ll(T,X)} = / / Y[p 2 (T i+1 - T u X i+l - Xi)e 2X ° l Xo >odX TT dTidX { 
JA> n (T) J K "+i f = \ 

Rescaling T { = TTi and X { = VTXi, we see that, 

E[ll{T,X)\ ~c„rt 

with C n ~ C\/re!. This shows that as T \ 0, 

E[(Z(T,X) -I (T,X) -h(T,X)) 2 } < CT. 

Note that Io(T, X) is just the solution of the heat equation with initial data Z(0, X). Let J£"(0) be 
the ex-field generated by Z(0,X) = exp{-B(X)}l x >o, X>0. Given J^(0), h(T,X) is Gaussian, 
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with mean zero and covariance 

Cov(h(T,X),h(T,Y)\3?(0)) 

= / / / p(T-T u X- X 1 )p(T 1 ,X 1 - X )e- B ^dX 

JO J-od \Jo J 

x (J"p(T - Ti,y - X 1 )p(T l ,X 1 - X' )e- B ^dX^j dX x dT x 

/*00 f'OQ 

= / / r- 1 /2*(T- 1 /2x,T- 1 /2y ) r-i/2 Xo ,r- 1 /2^).e- B ( Xo )- B ( x o)dx dxi 

Jo io 

Hence if we let 

Z init (T,X) = T-^h^T^X) 
it has the desired properties. □ 
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